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An Agricultural Chemical Transport Model 



By M. H. Frere, 1 C. A. Onstad, 2 and H. N. Holtan 3 



SUMMARY 



The use of chemicals is an important feature of 
modern agriculture, but the movement of the chemical 
from its point of application constitutes a potential 
hazard. The evaluation of this hazard is very difficult 
because the agricultural watershed contains many com- 
plex processes that dynamically respond to weather 
elements. To aid in evaluating mis hazard, a mathe- 
matical model has been developed and a computer 
program written that describes the movement of chemi- 



cals in and across an agricultural watershed. For each 
storm in series the objective is to use the model to 
predict the concentration of a chemical m the runoff 
water, the total amount carried by the runoff water and 
sediment, and the location and concentration of the 
chemical remaining on the watershed. This report 
describes the rationale used to develop the model as well 
as a couple of tests of the model mat illustrate and 
evaluate Its use. 



INTRODUCTION 



Maintaining a high quality of water resources, con- 
serving land resources, and increasing agricultural ef- 
ficiency are vital ingredients to national sustenance. The 
movement of chemicals and soil from agricultural areas 
can endanger these ingredients, 

The Federal Water Pollution Control Act Amend- 
ments of 1972 stress several requirements concerning 
agriculture. These include guidelines for Identifying and 
evaluating the nature and extent of pollution from 
agriculture lands, and the processes, procedures, and 
methods for controlling it. The Administrator of EPA, in 
cooperation with the Secretary of Agriculture and other 
agencies, is charged with developing the necessary 
research to provide the needed information. 



'Soil Scientist, Water Quality Management Laboratoiy, 
Durant,Okla. 74701. 
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Ignorance of essential factors in a complex system 
such as agricultural watersheds can lead to excessive or 
inadequate controls. The result can be economic or 
environmental damage requiring years to correct. The 
factors in a watershed system are interrelated, dynamic, 
and influenced by uncontrollable weather elements. Opti- 
mum management of an agricultural area requires 
predictions of the behavior of the watershed to changes 
in farm management practices. Quantification of the 
potential for any amount of chemical applied to any 
part of a watershed moving to the streams and ground 
water requires a comprehensive sequential accounting in 
time and space for the entire watershed. 

One of the better means for obtaining reasonable 
predictions of complex systems is to perform the 
accounting in terms of a mathematical model. The 
operation of the computational model on modern 
computers permits evaluation of numerous agricultural 
management alternatives and various climatic sequences. 
In addition, testing and development of a model serves 
as a research toot to integrate multidisciplinary projects 
and to reveal gaps in information and understanding. 



Hie essentials for computing movement of chemicals 
applied to agriculture require a distributed system of 
watershed hydrology, erosion, sedimentation, and chem- 
ical behavior. To establish a continuum for discussion 
and improvement, the USDAHL-74 Model was used for 
watershed hydrology, the Universal Soil Loss Equation 



(USLE) was adapted to compute the erosion, and 
sedimentation on physiographic areas in the watershed. 
A method was developed to compute the interactions of 
agricultural chemicals with the soil-water-plant system 
in the physiographic areas. 



OBJECTIVES 



The objectives for this model are to predict the 
concentration and amount of the chemical in the runoff 

water and on the sediment at the watershed outlet, and 
to predict the location and concentration of chemicals 
that are leached and moved spatially through the soil of 
the watershed, The model includes such management 
options as the time, the rate, and the type of chemical 

applied changes m crop pattern, and tillage practices. To 
facJitate application topographic, soil, and land-use 
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Considerable research will be required to achieve the 
objectives. This general framework constituting initial 
effoits divides the model into three submodels: A 
chemical, an erosional, and a hydrological submodel. 
The erosion submodel receives data fiom the hydrology 
submodel while the chemical submodel receives data 
from bom the hydrological and erosional submodels. 
This separation facilitates the interchange of other 
chemical, erosional or hydrological Submodels and modi- 
f ti nS - b * Ctives of the illitial effort "ere to achieve 
^ Perating V6rsion f AC for certain limited 
C nditi0nS - The model st then ^ tested, modified, 
Venf ' ed ' 3nd extended to <"* ^cultural areas before 
itcanbeco ^red a finished product. 



INITIAL MODEL BOUNDARIES 
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development because i, is one of the major areas of 
heavy ferU.zer and pesticide use. The initial version was 
jot designed to be apphcab,e to the important agricul- 
tunl watersheds located in flat coastal o, delta regions 
or containing irrigated or forested land. 

Ita transport scale is divided into t wo phase, n. 
flrst , cherracal application, depadation, redistribution 
by cu twaaon, and nitrate .runeraHxation and uptake by 
(he plants. The second phase is the transportation of 
chemicals by water and sediment movement. 

The operation of the model through the first phase 
following appDcation of a nonexistent chemical is 
adequate to descnbe the potential acute toxicity prob- 
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herb.cides do affect crop ve7 an d tat 
hydrologic and erosion ' 
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fact, the characteristics of these problems 
different that an entirely new checrf 
requ.red 

The erosion submodel is limited to rill and 
components of sheet erosion How 1" 
characterist.es of rill ' 

not considered. Local 
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a watershed, but differences m credibility and texture of 
tile deposited material aie not considered for subsequent 
storms. Sediment from areas such as gullies, roads, and 
construction sites is not accounted for, although such 
sediment may be significant in downstream chemical 
transport. The interactions of fertilizers, soil condi- 
tioners, and other cultural activities with the soil 
properties known to affect erosion are not considered, 
nor is the time rate-of-change of properties consideied. 
The hydrology submodel computes average values for 
evapotran 1 piration, infiltration, and runoff for each zone 



in the watershed. Therefore, chemical considerations 
limit a single crop for each zone because average values 
for dissimilar crops will not provide adequate informa- 
tion for chemical movement. Only two layers per zone 
are physically defined and provide useable data. Two 
other layers are empirically defined as storage from 
recession flows. The subsurface movement of water 
laterally from each zone is routed across the alluvium or 
directly into the channel. Snowmelt is computed in the 
model for all precipitation increments identified as snow 
in the input data. 



THE GENERAL FRAMEWORK 



As previously pointed out, ACTMO is divided into 
three sections to provide a means where alternate 
chemical, erosional, or hydrological submodels could be 
easily substituted. Figure 1 summarizes the general 
framework of the model in a flow diagram while table 1 
describes the information transmitted in the input and 
output files of ACTMO. 

The USDA Hydrograph Lab model was adopted as 
me hydrologic submodel because it was operational, was 
calculating most of the variables considered necessary 
for a transport model, and utilized parameters readily 
available from soil suivcy data. The hydrology submodel 
was developed from data on a variety of watersheds 
across the country and, therefore, is potentially broad 
based. The eiosional and chemical submodels were 
developed specifically for this model, combining many 
previously published concepts. 




The processes involved m chemical transport indicate 
tile following factors are important. The soil moisture 
and temperature influence the breakdown of many 
chemicals and the mineralization of nitrogen, but nitro- 
gen uptake is expected to be correlated with transpira- 
tion. The leaching efficiency of water tends to be lower 
as gravitational water percolates through the soil. The 
amount of runoff is obviously an important variable as is 
the relative amount of lateral flow for mobile chemicals. 
Therefore, the hydrology submodel supplies the chemi- 
cal submodel the following variables; The temperature, 
and for each layer in a zone the evapotranspiration and 

Table I. -Description of the information transmitted in 
the input-output files of ACTMO 



File 



Information 



Figure l.-Flow diagram for the input-output flies In ACTMO. 



Weather , . . Rainfall, temperature, pan evaporation 

Watershed- Zone sizes and cascading, soil hydrologic proper- 
ties, crop pattern and characteristics, tillage 
operations. 

Storm .... Length, slope, and cascading of zones, runoff, 
infiltration, und rainfall intensity. 

Soil Soil texture, 

Erosion .... Compartment areas, enrichment coefficient, soil 
texture, infiltration, runoff, erosion, and deposi- 
tion. 

Hydrological Cultivations, soil properties, time-averaged tem- 
perature, soil moisture, and ET, lateral and 
vertical flow, soil water when the storm starts. 

Chemical , . Date, amount, and type of chemical application; 
treated area; soil chemical interactions; and 
critical levels. 

Output of Chemical lost in runoff and erosion. 
Chemical Concentrations in runoff. ' 
Behavior. , . Chemical left In each compartment. 

Location and concentration of peak in each 

compartment. 



soil moisture averages between storms, the infiltration 
before runoff starts, infiltration during runoff, volume 
of runoff, the lateral and vertical flow rates for each 
zone, and the soil moisture in each layer of each zone at 
the start of the storm. 

Chemicals adsorbed to the soil are subject to erosion 
as a transport mechanism. The chemical concentration in 
the soil varies with depth because of the application 
method and the amount of leaching. Therefore, both rill 
and interrill or surface erosion are calculated. Finally, a 
decrease in the transport capacity of the runoff can 
cause deposition of the sediment. Chemicals are usually 
preferentially adsorbed by the finer sediment fractions 
which are the last to be deposited and, therefore, the 



composition of the sediment which is enriched with the 
fine material is calculated. In summary, the erosion 
submodel supplies the chemical submodel with estimates 
of interrill erosion, rill erosion, area and depth of rills, 
amount of deposition, and the enrichment of the eroded 
sediment relative to the clay fractions of the original 
soil. 

The Universal Soil Loss Equation was modified to 
compute the erosion and deposition information for the 
chemical submodel. The rainfall factor in this equation 
was adjusted from an annual average to a storm basis by 
the rainfall intensity, the amount of runoff, and the 
peak rate of runoff variables supplied by the hydrology 
submodel. 



THE WATERSHED TRANSFORMATION 



On a multicropped watershed one must anticipate 
that only part of the watershed might be treated with a 
chemical. Some means must be used to transform the 
watershed defined by physical characteristics, such as 
crop fields, into computational units. A zone, the 
computational unit for the hydrology submodel, is 
constructed by grouping fields of the same crop, or 



major soil areas, or land-capability classes or whatever 
physical feature is recognized as important. The zones 
typify the elevation sequence of uplands, hillsides, and 
bottomlands. The percent of each zone contributing 
surface runoff to lower zones is estimated by observa- 
tions from the soils and topographic maps. The average 
downslope dimensions of each zone is the length of 
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Figure 2. -A three-dimensional transform of the watershed into computational units. 

4 



overland flow in the zone. The results are illustrated in 
figure 2 where foui zones with different slopes are each 
divided into four compaitrnents creating seven tubes. 

Each tube is the area between adjacent streamlines 
from flow net theory. The flow of water and sediment is 
assumed to occur along the line of steepest descent, 
representing flow over an i&otropic surface. This provides 
a situation analogous to the natural sequence of flow 
while maintaining the distributed effect. The size of the 
incremental areas within a zone or compartment is 
calculated from the average zone length and the fraction 
of the zone area that contubutes directly to the channel. 
Each compartment is intended to be a homogeneous 
unit for computing infiltration, evapotranspiration, over- 
land flow, and net erosion or deposition. Soil and water 



lost from one compartment are considered to be 
additions to the top of the adjacent downslope compart- 
ment wheie similar calculations are made. This process is 
continued down to the watershed mam channel. 

Figure 2 also illustrates the vertical dimension that is 
considered. The hydrology submodel treats the first 
layer as the "A" hoiizon 01 plow layer and the second 
layer as the remainder of the potential loot zone. The 
other two layers are defined empirically as storage from 
recession outflows. The chemical submodel uses the 
same layers to compute the chemical interactions with 
the soil. Although coefficients of adsorption, dispersion, 
and breakdown could be handled as continuous func- 
tions of depth, the degree of lefinement is limited by the 
degree of hydrological definition of the soil. 



SUBMODEL DESCRIPTION 



This section provides a detailed description of each 
submodel. Such a description is not essential for an 
overview, but a model, like a chain, is no better than its 
weakest link. In this section all of the iclations and 
assumptions are laid bare for examination. 

The description of each submodel is essentially an 
independent paper. The section on the hydrology 
submodel is a recent revision of "USDAHL-70 Model of 
Watershed Hydrology" (20). 4 The erosion-deposition 
submodel was developed by C. A, Onstad while M. H. 
Frere developed the chemical and nitrate submodels. 

The Hydrology Submodel for ACTMO 

The USDAHL-74 Revised Model of Watershed Hy- 
drology 5 

Hydrologic Grouping of Soils and Land Use 

Soils on each watershed are grouped by land capa- 
bility classes, as illustrated in figure 3, to form hy- 
drologic response zones for computing infiltration, 
evapotranspiration, and overland flow. Zones indicated 
for each watershed typify the elevation sequence of 
uplands, hillsides, and bottom lands in these areas. Zones 
are always numbered in a downslope order because 
computations assume that some of the runoff will 
cascade over successive zones. For example, on the 
Coshocton watershed an estimated 60 percent of Zone 1 
runoff flows directly onto the alluviums of Zone 3, and 
only 40 percent cascades to soils of Zone 2. 



^Italic numbers in parentheses refer to Literature Cited, 
p. 37. 

5 USDA-ARS Plant Physiology Institute Report No. 4 by 
H. N, Holtan, G. J. Stilner, Wm. H. Henson, and N. C. Lopez, 
1974. 



The watershed description data published by the 
Agricultural Research Service field stations (48} give 
percentages of area of major soils and ofland-capabihty 
classes. Moisture-tension data (table 2) published by the 
USDA Hydrograph Laboratory (22) or the Soil Survey 
Investigations Reports of the Soil Conservation Service 6 
provide soil profile descriptions for computing weighted 
averages of hydrologic capacities (8) of soils in each 
zone. Specific gravity was estimated at 2.65 for all soils 
in computing total porosities. Moisture at 0.3 bar 
tension was used for field capacity in dividing total 
retention capacities (S) into moisture freely drained by 
gravity (G) and moisture drained by vegetation (AWC). 
The lower limit of AWC is assumed to be the moisture at 
15 bais tension. Table 3, by C. B. England, is a summary 
of a great many samples giving S, G, and AWC in percent 
volume of textural horizons. Percentages are applied to 
depths of "A" horizons or to plow depths in nonlayered 
soils, and to depths of maximum root penetration as 
might be limited by hardpans, gravelly layers, or other 
textures unsuited for root growth. The input format of 
the model accepts G and AWC in percentages by volume 
and applies them to soil depths to compute volumes. 

Evapotranspiration 

Evapotranspiration potentials are estimated by coef- 
ficients applied to published pan-evaporation data. The 
method is a combination of techniques developed by 
Mustonen and McGuiness (52), Jensen (25), and 
Pruitt (38) adapted for use with cardinal temperatures 
presented in table 4 in estimating plant growth. Figure 4 
is a plotting of XGI, cardinal temperature functions 



6 Soil Conservation Service. Soil Survey Investigations Report 
Series by States, U S. Dept. Agr., Soil Conserv, Serv. 



Table 2 -Profile description and physical properties- Cecil Gritty Sandy Loam 

Location Blacksburg, Virginia, Brunswick County, Rocky Run Branch Watershed, located in pasture 75' SW of 

fence, 400 yds. W of county road #656. 
Vegetation and land use: Pasture. 
Topography. Gently sloping. 
Drainage Well-drained 
Parent material' Granite gneiss and schist. 
Described and sampled by. K. Fussell, J. Williams, and J. B. Burford. 



Horizon 



Description 



Available moisture 



per inch of 

soil 



In horizon 



Ap, 



B22t.. . 



B231.. .. 



B31 



B32 



to 6 inches. Dark yellowish brown (10YR 4/4) gritty sandy loam, very friable; weak 
medium granular structure; many fine roots, medium acid; few small peds of red (2.5 YR 
4/6) brought up from lower horizons by plowing; abrupt wavy boundary. 

6 to 15 inches. Red (2.5YR 4/6) clay with few fine mottles of yellowish red (5YR 5/8), 
some ped faces are coated with yellowish red (SYR 4/8) clay, which is apparently result 
of root and worm action; firm; moderate medium subangular blocky structure, thin 
continuous clay films; few mica flakes; common fine roots; occasional quartz rock less 
than 1/2" in diameter; medium acid; clear wavy boundary. 

15 to 31 inches. Red (2 SYR 4/6) clay loam with common fine faint mottles of yellowish 
red (SYR 5/8) and common fine distinct mottles of yellowish brown (10YR 5/8); moderate 
medium and fine subangulai blocky structure; friable; thin clay film on most ped faces; 
common mica flakes, 10 percent grit by volume, few fine roots; strongly acid; clear wavy 
boundary. 

31 to 42 inches. Red (2. SYR 4/6} clay loam with many medium faint mottles of yellowish 
red (5YR 5/8) and few medium distinct mottles of strong brown (7. SYR 5/8) and reddish 
yellow (7.5YR 6/8), fnabte; moderate to weak medium and fine blocky and subangular 
blocky structure; patchy clay film; many mica flakes; 10 percent grit by volume; strongly 

acid , clear smooth boundary. 

42 to 56 inches. Mottled colors of red (2. SYR 4/8), yellow red (SYR 5/8), reddish yellow 
(7. SYR 6/8), very pale brown (10YR 7/4) highly weathered granite gneiss and schist with 
few pockets of light clay loam , friable, weak fine and medium subangular blocky structure; 
many mica flakes; slimy and greasy to the feel, 10 percent grit by volume; few patchy clay 
films; ver> strongly acid;gradual wavy boundary. 

56 to 71 inches plus Mottled colors of yellowish red (SYR 5/8), red (2 SYR 4/6), reddish 
yellow (5YR 6/8), and very pale brown (1 OYR 7/4) highly weathered granite gneiss and 
schist; silt loam with 10 percent grit by volume. 



Inches 
0.10' 



.07 



.06 



.13 



.18 



.18 



Inches 

0.60 



.63 



.96 



1.43 



2.52 



Horizon 


Depth 


Bulk 
density 


Moisture as percent of dry weight at respective atmospheres of 

tension 


01 


0.3 


0.6 


3 


15 


An 


Inches G Ice. Percent 
0-6 1.72 15.24 
6-15 1.33 26.90 
15-31 1.32 32.02 
31-42 1.41 29.89 
42-56 1.29 35.15 
56-71+ 1.29 36.36 


Percent 
12.35 
23.76 
29.26 
28.02 
26.66 
24.46 


Percent 

10.30 
22.22 
28.67 
25.87 
19.52 
26.29 


Percent 
9.21 
21.53 

27.27 
21.68 
16.49 

15.78 


Percen t 
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18.28 
24.58 
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13.02 
10.69 
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Figure 3.-Hydrologic zones of land-capability classes, 



Table 3.~Hydrohgic capacities of soil texture classes 

[Adapted from "Land Capability: A Hydrologic 
Response Unit In Agricultural Watersheds," by C. B. 
England, Res. Soil Scientist, USDA Hydrograph 
Laboratory, Beltsvllle, Md. (USDA ARS 41-172, 
Sept. 1970)] 



Texture class S 1 G 2 AWC 3 AWC/G 



Percent Percent Percent 


Coarse sand 


24.4 


17.7 


6.7 


0.38 


Coarse sandy loam 


24.5 


15,8 


8.7 


0.55 


Sand 


32.3 


19.0 


13.3 


0,70 


Loamy sand 


37.0 


26.9 


10.1 


0,38 


Loamy fine sand 


32.6 


27.2 


5.4 


0.20 


Sandy loam ~ . . . . 


30.9 


18.6 


12,3 


0,66 


Fine sandy loam , 


36.6 


23.5 


13,1 


0.56 


Very fine sandy loam 


32.7 


21.0 


11.7 


0.56 


Loam 


30.0 


14.4 


15.6 


1.08 


Silt loam 


31.3 


11.4 


19.9 


1.74 


Sandy clay loam 


25.3 


13.4 


11.9 


0,89 


Clay loam 


25.7 


13.0 


12.7 


0.98 


Stlty clay loam 


23.3 


84 


14.9 


1.77 


Sandy clay 


19.4 


11.6 


7,8 


0.67 


Silty clay 


21.4 


9.1 


12.3 


1.34 


Clay 


18.8 


7.3 


11.5 


1.58 



S = Total porosity - 15 bar moisture percent. 
G= Total porosity - 0.3 bar moisture percent. 
AWC = S minus G. 



(current temperature - lower cardinal temperature)/ 
(optimum cardinal temperature -lower cardinal tempera- 
ture) for corn and alfalfa. Cardinal temperatures are the 
optimum and lower or minimum temperatures for crop 
grow til. When current temperature exceeds the optimum 
cardinal limit, the plant is assumed to suffer and the XGI 
function is set less than 1 .0 by one-half of the amount of 
exceedence. If current temperature is equal to or less 
than the lower cardinal limit, the XGI is set equal to a 
fixed minimum of 0.1 . 

The growth index for a given crop is then computed 
from the XGI curve for that crop, as illustrated in figures 
5 and 6 for corn and alfalfa respectively. The crop 
growth index, GI, follows the XGI if no cultural 
practices intervene. However, cultural practices that 
reduce the foliage also reduce evapotranspiration; hence, 
GI for the crop is reduced following plowing, planting, 
cultivating, or harvesting if vegetation is abolished or 
significantly reduced by the operation. (This excludes 
minimum tillage practices.) After each practice, the crop 
is assumed to recover in a specified number of weeks at 
which time the GI curve rejoins the XGI curve. In the 
case of turnplowing, the soil is inverted and left in an 
aerated condition; hence, the GI value is set at 1.0 and 
assumed to settle back to a bare fallow condition of low 
ET in two weeks where it remains a constant at 0.1 until 
planting of a crop. 



Cardinal temperatures can be obtained from the midsummer, the associated temperatures indicate the 
literature for most crops. Table 4 is a compilation of upper limit. Winter wheat would be expected to have a 

lower limit near to freezing. The evapotranspiration 
potential can then be estimated for any given day as 



upper and lower cardinal temperatures for some of the 
more common crops together with literature citations. 
Generally, the local reputation of a crop for seasonal 
growth can be drawn upon to approximate the cardinal 
limits; e.g., if a grass is known to go dormant at 



k- E [(S-SA)/AWC] 



[1] 



Table 4. -Cardinal temperatures for growth of some common crops, degrees Fahrenheit 



Crop 


Cardinal temperatures 






Minimum Optimum Maximum 




Oats, rye, 
wheat, barley. , 

Sorghum, 


32-41 77-87.8 87.8-98.6 

59-64.4 87.8-98.6 111.2-122 

50 
40 
55 

40 
50 86-95 113 

49 93 115 
35 70 82 

33-34 

46-50 
55-57 

no 

50 60-95 105 
60 
80-84 
100 
60 80 100 
60 


Parker/N. W. 1946. Environment factors and their control in 
plant environments. Soil Sci. 62:109-119. 

Chang, J. H, 1968. Climate and Agriculture: An Ecological 
Survey. Aldme Pub!. Co., pp. 75-76 

Wiggins, S. C. 1936. The effect of seasonal temperature. . . 
Agron. J. 48. 21-25. 

Meyer, B. S. and D, B. Anderson, 1939. pp. 589-591. D. van 
No strand Co., Inc., New York. 

Costing, H. J. 1953. The Study of Plant Communities. W. H. 
Freeman & Co., San Fran., pp. 127-128. 

Haberlandt, 1874, cited by Richards, S. J. et al 1952. Soil 
temperature and plant growth. In Soil Physical Conditions 
and PI,mt Growth. Agron. Mono. II: 360-363. 

Lehenbauer, P. A. 1914. Growth of maize. . .Physiol. Res 

1; 247-288. 

Wang, J. Y. 1963 Agricultural Meterology. Pacemaker Press, 
Milwaukee, pp. 104-105. 

Rhykerd, C. L., et A\. I960. Sorghum grows best at warm 
temperatures. Crops and Soils 12-24. 

Martin, J H. 1941. Climate and Sorghum, in USDA Yearbook 

of Agric , p. 343 

Doggett, H. 1970. Sorghum. Longmans, Green and Co., Ltd. 
London, pp. 188-189 

Garner, W W. 1946. The Production of Tobacco. Blakiston Co 
Phila., pp 376-377. '' 

Weihing, R. N 1963. Growth of ryegrass as influenced by 
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Figure 4,-Rclativc potential growth (XGI) curves for corn and alfalfa computed from upper and lower cardinal 

temperatures (TU and TL), 
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Figiiro 5. -Relative growth (GI) curve for corn as modified by tillage, 
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I igurc 6 Relative growth (GI) curve for alfalfa as modified by cutting 



H"i evapotranspiration potential in inches per 

day , 

GI is growth index of crop m percent of maturity, 
k is ratio ul dl to pan evaporation, usually 
1 0-1.2 fur short grasses, 1.2-1 6 for crops up to 
shouldt-'r-IiugJit, 3 6-2 lur forest, 
I-p b jun evaporation m inches per day. 

Si) roul porosity, 
SA LI available porosity, 

AWf is porosity dramable only hy evapotranspira- 
unu, 

and \ is set equal to 0.5 for sands and 1.0 for finer 

textures 

Input to (hi; model now includes 52 weekly averages 
ot air temperature in plate of tin,' Gl curves required for 
ejt'li crop in Model-70 The use of temperature is 
designed to individualize plant growth estimates for each 
uMr, 



Dr Wadleigh, m mathematical analyses of plant 
growth data (unpublished), demonstrated that free water 
(above field capacity) dampens and inhibits plant 
growth, and subsequently ET, by exclusion of oxygen 
needed for plant growth. As illustrated schematically in 
figure 7, growth and, consequently, ET decreases from 
optimum at field capacity to zero at soil satuuition. 
Thus ET increases from wilting point to an optimum a t 
field capacity and diminishes to zero at saturation. 
However, free-water evaporation increases from zero at 
field capacity to a value equal to pan evaporation at soil 
saturation. Tins separation of evaporation and transpira- 
tion becomes important in chemical transport studies. 

Infiltration 

Infiltration capacity was expressed by Holtan(/5) as 
an exhaustion phenomenon convergent upon some 
constant rate" 



f=a-S 



14 



[2] 
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where 

f is infiltration capacity in inches per hour; 
a is the infiltration capacity m inches per hour per 
(inch) li4 of available storage (index of surface- 
connected porosity); 

S a is available storage in the surface layer (the "A" 
horizon in agricultural soils) in inches water 
equivalent; 

and f c is the constant rate of infiltration after pro- 
longed wetting (associated with capillary flow or 
with an impeding stratum) in inches per hour. 
The GI of equation [1] is also used as a seasonal 
factor on the vegetation parameter "a" in the infiltration 
equation [3] : 



f=GI- a- S 



1.4 



[3] 



Infiltration and rainfall excess are computed for each 
soil zone by comparing observed rainfall to the infiltra- 
tion capacities computed by equation [3]. In solutions 
of equation [3] , the volume of infiltration cannot 
exceed the volume of water available at the soil surface 
(including detentions of overland flow), nor can it 
exceed the product of the time increment and the 
average rate, (f t + f 2 ) * t/2. Infiltration capacity de- 
pends on the amount of empty space with S a diminished 
by infiltrated water but recovered by drainage or 
evaporation to the limit of G, and by ET to the 
additional limit of AWC. 

Values of "a" are shown in table 5 for various land 
uses in agriculture. These were evaluated at plant 
maturity as the percentage of the ground-surface area 
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MATRIX SUCTION IN BARS 



Figure 7. -Relative growth as modified by free water (matrix 
suction less than .25 bus) and moisture deficit (matrix suc- 
tion greater than .25 bars), (After Wadleigh, unpublished.) 



Table 5. - Tentative estimates of vegetative parameter 
"a" in infiltration equation 



Land use or cover 


Basal area rating 1 


Poor 
condition 


Good 
condition 



Fallow^ 


0.10 


0.30 




.10 


.20 




.20 


.30 




.20 


.40 




.40 


.60 




.20 


.40 




.40 


.60 




.80 


1.00 




.80 


1.00 









Adjustments needed for "weeds" and "grazing." 
For fallow land only, poor condition means "After row 
crop," and good condition means "After sod." 



occupied by plant stems or root crowns. Stems and root 
crowns are assumed to reflect the percentage of porosity 
in the top layer that is surface-connected by mature 
plant roots to form conduits for air or water. If surface 
continuity is achieved by incorporation or by me- 
chanical means, the "a" value should be adjusted to 
represent it. 

Estimates of f c were obtained from the SCS Hand- 



book (45) wherein major soils of the United States are 
grouped according to their rate of water intake after 
prolonged wetting. Soils are grouped in four hydrologic 
classes; A, B, C, and D. Musgrave (30) gives the as- 
sociated rates of f c in inches per hour: A = 0.45 to 0.30; 
B = 0.30 to 0.15; C = 0.i5 to 0.05; and D = 0.05 to 
0.00. The texture and density of the impeding layer 
give a clue to the selection of f c within a group. If the 
impeding strata approach clay, f c is near the lower 
limit of its group; for sand f c would be near the top. 
Loams would indicate the midpoint f c . 

Depressions on the land surface act as reservoirs 
holding a volume, V d , of water until it is dissipated by 
infiltration. Although generally quite small in depth 
(about 0.05-inch for agricultural crops), Vd is drawn 
upon by infiltration and subsequently refilled by rainfall 
for each pause in rainfall intensity. Also, certain agricul- 
tural practices such as level terraces, contoured furrows, 
and listed crops may be designed to give Vd depth of 
several inches. Musgrave and Holtan (32) give concepts 
and some estimates of depression storage in agriculture, 
but generally such, estimates must be based upon 
personal observations. 
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Certain soils, such as the montmorillonitic clays, form 
deep cracks on drying. Cracking is estimated for a given 
horizon from the ratio of bulk density at field capacity 
to bulk density when air dry by the equation: 



% Cracks = 100 



BDD 





BDD 



[4] 



where 



BDWis bulk density at field capacity in percent 

volume; 
and BDD is bulk density at air dry in percent volume. 

Tliis equation was developed by considering a cube of 
soil and by deleting vertical shrinkage. Equation [4] is 
applied to layer thicknesses to compute maximum 
cracking volumes. The volume of cracks at any given 
time is computed as a linear function of soil moisture 
present, limited to the plant available range, AWC. 
Cracks are therefore at maximum within the root depth 
at wilting point and disappear at field capacity. 

Coefficients for Routing Flows 

Rainfall in excess of infiltration is routed across each 
soil zone of figure 3 and cascaded, subject to further 
infiltration, across designated subsequent soil zones 
enroute to the channel. Overland flow is computed by 
an adaptation of the continuity equation: 



and 



where 



q = (ova)D 



1.67 



15] 



16] 



Peis the volume of rainfail per At in excess of 

infiltration and depression storage; 
Q is volume of outflow per At; 
QO is overland flow in inches per hour; 
ova is a coefficient dependent upon roughness, length, 

and degree of slope; 
D is the average depth of flow in inches; 
and t is time in hours. 

The coefficient "ova" in equation [6] is determined 
from runoff recessions on rectangular plots as described 
by Musgrave and Holtan. It is the rate of flow produced 
by an average depth of 1.0 inch of water on the plot 
surface. In the flow data presented by Musgrave and 



Holtan, "ova" can be estimated from vegetative density 
"a", percent slope "SI", and feet length "Le" of flow 
path by the following equation: 

ova = (150 - 125a)-(12/ VO) ' VsT/Le [7] 

For laminar flow, "ova" would vary linearly as "SI" in 
equation [7] and the exponent, 1.67, must be changed 
inequation [6] to 3.0. 

Channel flows and subsurface return flows are routed 
by simultaneous solutions of the continuity equation 
and a storage function. Storage coefficients are obtained 
by integration of the flow recession curve for a given 
watershed. Techniques are under study in the USDA 
Hydrograph Laboratory for the approximation of flow 
recession curves by linear segments to obtain the 
advantages of superposition in routing water flows from 
different areas of the watershed. Flow from each unit is 
routed separately through watershed storage, and all are 
summed to obtain watershed outflow. 

Postrainfall recession curves are plotted on semi- 
loganthmic paper in figure 8 for Little Mill Creek 
Experimental Watershed near Coshocton, Ohio. Evapo- 
transplration must be minimized if we are to obtain a 
true storage-flow relationship. The winter (March) re- 
cession shows the least influence by ET and is used for 
analyses. The equation of a recession curve segment is: 



-t/m 



[8] 



where: 



q is rate of flow at start of a period in inches per 

hour, 
q t is rate of flow one time-increment later in inches 

per hour, 
m, the absolute value of t/Aln(q), is a constant for 

each straight-line segment of the recession curve 

on semilogarithmic papers; 
e is logarithmic base; 
and t is the time-increment in hours. 

By integration of equation [8] the storage increment, 
AS, within a linear segment of the recession curve is: 



AS = m Aq 



[9] 



Values of "m" derived for each linear segment of the 
semilogarithmic plotting are assumed to represent suc- 
cessive flow regimes (34 J, starting with m c for channel 
flow and proceeding through a series-"m i , m 3 , m 3 , and 
m 4 "-for successively deeper or more devious regimes of 
subsurface flow. Some watersheds, such as the Hastings 
W-l 1 , have no return flow; hence, only m c is defined. 
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Figure 8. -Seasonal recession flows, W-97, Coshocton, Ohio. 



At Coshocton, flow regimes must be separated. The 
usual procedure has been to extend segments backward 
in time for subtraction from earlier regimes (5). In our 
experience, and that of Kirlandaiswamy and Seethor- 
man (27), the adjusted points never plotted as a straight 
line, and the resulting recession was too rapid. The 
authors prefer to use the same procedure generally 
used for separation overlapping hydrographs of surface 
flows, that is, extending the recession of the first storm 
forward in time for subtraction from subsequent flows. 

Figure 9 illustrates results of techniques suggested for 
separating flow regimes for input to this model. Each 
linear segment ,is extended toward zero and subtracted 
from subsequent flows to get the rising hydrograph of 
the succeeding regime. Now there is no question about 
maintaining logarithmic linearity within flow regimes. 
Also, the rising hydrograph fits better with the concept 
that return flows are less when flows in stream channels 
are at higher stages. The total volume of the separated 
hydrograph can be" computed for successive subsurface 
regimes by the equation: 



[10] 



where: 

S is total storage volume of the difference hydro- 
graph In inches; 
q b is rate at intersect with the previous segment of 

recession; 

m is absolute value of t/ln(q /q t ); 
L is flow regime number; 

and subscript L-l indicates "m" of the previous 
segment. 

Input to the model would then be "m c , qi, raj, q^> 
ma, qa, m 3 , q 4 , and m 4 > ' for those watersheds having 
four discernable regimes of subsurface flow, The value of 
'W is readily determined as the number of hours 
required for the recession segment to cross one log cycle 
divided by 2.3, the natural logarithm of 10.0. Maximum 
free-water storage is assumed equal to the "G" value 
from the soil survey for regimes 1 and 2, but it is the 
volume under the hydrograph computed by equation 
[10] for the remaining flow regimes as in figures 8 and 
9. 

M c , as determined from the watershed flow recessions, 
Includes detentions of overland flow. Since overland flow 
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SEPARATION OF FLOW REGIMES 
W-97, Coshocton, Ohio 
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Figure 9. -Separation of flow regimes, \V-97, Coshocton, Ohio. 



is routed in equation [6], the depth of overland flow 
must be subtracted from the AS of equation [9] . At the 
peak, outflow is essentially equal to inflow and channel 
storage would be. 



= m c . q~(q/ova)" 6 



[11] 



We can get an approximation of the corrected m c 
(cm c in the program) by solving for m in equation [9] at 
one unit of flow rate. 



cnv - 



ova 



[12] 



Equation [12] is an approximation only and typifies 
an area of needed research. 

Concern For the transport of agricultural chemicals 
has placed emphasis on the path of both surface and 
subsurface flows from relatively small upland areas well 
above the elevation of emerging subsurface return flow. 
As a result, values of m v to m 4 are obtained from 
downstream gaging sites of an encompassing watershed 
together with an indication of the number of regimes 
(including m c ) which occur above the weir. Output then 
includes "onsite" return flow in "runoff and lists 
"offsite" flows for comparison with regional informa- 
tion. 



As the size of watershed diminishes, the importance 
of detail increases. Therefore, mi and ma are subscripted 
by zone and Devaluated from soils information. As- 
sumptions were: 

(1) Outflow is equal to permeability width/length 
gradient; 

(2) permeability varies as the ratio of pore diameters 
squared or as the square of the ratio of large pore 
volumes; 

(3) and, since f c is permeability of layer 2> the 
permeability of layer 1 is f c (Gi/G a )*. 

Then, since storage, (FREE WATER) is equal to 
DEPTH/G and since m ~ FREE/G, we can derive mz(l) 
and mz(2) for each zone by equation [13] , 



mz(L) = depth (L) G 2 3 Le 2 (SI 2 + 1)' 
(f G(L) zone area SI) 



[13] 



If the downstream flow record indicates that the 
hydrograph recession passes abruptly from the character- 
istic surface flow hydrograph to the low, sustained 
baseflow without transitional slope breaks, mi and m 2 
can be entered as zeroes. The program will not compute 
flows in a regime whose m value is zero. 
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Equation [14] is a combination of the continuity 
equation and equation [9] used in computing outflow 
from each regime: 



2AI 



2m -At 
2m +At 



[14] 



Equation [9] , being linear, forces the routed hydrograph 
from equation [14] to peak at its crossing with the 
inflow hydrograph, whereas channel flow responds non- 
linearly to changes in inflow rate. Therefore, non- 
linearity was induced on the rising side of the 
streamflow routing by the techniques of Hoi tan and 
Overton (21) for successive routings through one-half of 
indicated storage using equation [15] : 



, 2AI 




in -At 



m + At + qi m + At 



[15] 



In the second application of equation [15] the quantity 
(qi '+ Q2 ')At/2 is substituted for AI, The symbols are: 

I is inflow volume in inches; 
m is routing coefficient in hours; 
t is increment of time in hours; 
q is rate of outflow from storage in inches per hour; 
q' is theoretical rate of outflow from first half of 

storage jn inches per hour; 

and subscripts i and 2 indicate beginning and end of 
At. 

Hydrogeology 

In the USDA Hydrograph Laboratory's model, in- 
filtrated water is proportioned to ET, or free-water 
evaporation, and to downward seepage, or to lateral 
return flow, in each flow regime. Downward seepage and 
lateral flow are supplied by free water; hence, estimates 
of maximum seepage rate "C" and free water capacity 
"G", as well as the storage coefficient "m", are needed 
in each flow regime. 

Since subsurface flow regimes are considered se- 
quential, seepage, C, from a given regime is inflow to the 
next regime and must tre adequate to supply the sum of 
the maximum flow rate "q. " experienced in the next 
regime and the rate of deep seepage, gr: 



1 L+1 



[16] 



where: 



L is the regime number; 

and: gr is the maximum rate of ground water recharge 
in inches per hour. 



Ground water recharge from the ultimate return-flow 
regime is estimated on a regional basis Average annual 
rainfall, average annual ET, and average annual stream- 
flow yields in the vicinity can be used to derive an 
average annual ground water recharge. This can be 
converted to inches per hour as an estimate of gr. At 
Hastings, this value was estimated as equal to f c ; but the 
estimate for the Coshocton, Ohio, watershed was 0.0006 
inches per hour, and for the Texas and Florida water- 
sheds, gr was estimated as zero. 

It must be borne in mind that for small upland areas 
the deeper regimes may pass beneath the point of 
surface-flow measurement. If so, gr must be estimated 
equal to q max of the first unmeasured regime and the q 
and m values of the nonoccurring regimes must be 
omitted. 

Increments of downward seepage, "subput", to the 
next regime, are computed as a function of a free-water 
present: 



Subput L+1 = t C (G - SA)/G 



[17] 



The potential rate of outflow q p from a regime is 
derived by equation [14], but actual outflow from a 
regime is held equal to zero, as illustrated in figure 9, 
whenever summation outflow exceeds the value of 
intercept q b for that regime. When summation outflow 
is less than q b , a favorable gradient is assumed and 
outflow q L is computed by the following equation: 



{18] 



where 



q b is intercept with previous segment; 
2^ is sum of all existing flow rates; 
and q_ is potential regime outflow from equation 
[141. 

An exception to equation [18] is made for frozen 
ground. If the 3-week moving average air temperature is 
below freezing, the c of equation [17] is reduced for 
layer 1 by 0.02 Inch per hour for each average degree 
Fahrenheit below 32. Thus It is assumed that highly 
permeable soils become Impervious at lower tempera- 
tures. Since only f c is affected, soils take in water up to 
saturation during low temperatures but no outflow 
occurs from layer 1 until temperatures rise above 
freezing. 

The limits of free-water storage in the first two 
regimes are computed from soil survey data as the 
products of horizon thickness and percent freely drained 
porosity. The limit of free-water storage, G, associated 
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with subsequent flow regimes is computed from the 
storage-flow relationship: 



[19] 



There is no assurance that the flow recession analyzed 
represents a saturated A or B horizon. Hence, the 
G-v allies determined from soil survey data for these twc 
honzons take precedence in the model. Subsequently, 
the q-values and (lie G-values of regimes 3 and 4 are 
apportioned to each /.one of the watershed in proportion 
to the free-wa tcr limit "G" of the topsoil: 



G 3 (zonc) = 



lnia\3 



G, (/one)/G, [20] 



where GI is a weighted watershed average from soil 
survey. 

By (Ins means the storage indicated by the lower flow 
recession is dispersed over the watershed in proportion 
lo the depth of the topsoil horizon in each zone or 
hydrologic response unit. Routing coefficients, 1113 and 
ni4, remain unchanged for each regime, but the storage 
potentials, G n , and the maximum outflows, q raax , are 
now unique for each zone. 



[21] 



Evapotranspiration is drawn from both G and AWC 
when roots are present. Geologic research is needed to 
develop concepts of flow regimes associated with spe- 
cific depth in the profile. In the meantime, we are 
placing greatest icliance on soil surveys for storages of 
regimes 1 and 2 and depend on the recession curve 
analyses for storage values of subsequent regimes. Roots 
are arbitrarily designated as drawing water from the first 
two regimes, often creating a deficit in AWC to be 
satisfied before accruals to subsequent regimes can 
occur. Therefore, AWC| and AWC 2 are required input. 
AWC| is. obtained from the moisture tension data of the 
first layer. AWC 2 is estimated from the moisture tension 
data below the first layer to the depth of suitable habitat 
for roots as indicated in the soil profile description. 
Moisture deficits are computed only in the depth 
indicated for root habitat. Today, there is no basis for 
distributing any part of AWC to subsequent regimes 
except in a few places like Florida where records of wells 
over impervious rock at very shallow depths may permit 
some speculations. 

Evaporation from free water in the soil is treated as a 
function separate from evapotranspiration. II is ex- 
tracted from the sum of free water "G" in regimes 1 and 
2 as the product of percent total free water available and 
pan evaporation. 



A percentage of overland flow horn each one is 
designated to cascade across the subsequent soil zone, 
with the remainder, if any, allocated to the alluviums or 
directly to channel inflow Allocation is a matter of 
judgment based on the nearness of a zone to alluviums 
or upon the presence or absence of diversions, graded 
terraces, listed furrows, or other interceptors leading 
directly to a stream channel. Infiltration in excess of 
AWC in the root zone is routed through storage of 
subsequent regimes and accumulated as run-on to the 
alluvium or as direct inflow to the stream channel if so 
designated by input. The total of all inflows to the 
channel is routed, using the channel storage coefficient 
m c in equations [13] and [15] to obtain the outflow 
hydrograph for the watershed. 

The Erosion-Deposition Submodel 

"A Watershed Sediment Yield Model with 
Hydrological and Chemical Interfaces" 

This submodel predicts components necessary for the 
erosion-chemistry interface. It simulates the erosion- 
deposition continuum for watershed slopes on an area 
basis. This phase is subdivided to estimate the interrill 
and rill erosion components of the total sediment load. 
Particle size distributions of the eroded material are also 
estimated for erosion-chemistry relations. Each calcula- 
tion utilizes hydrologic input together with basic soil 
parameters for each compartment of the watershed 
transformation. 

Basic Equation 

The conceptual basis for the erosion and sedimenta- 
tion submodel is the Universal Soil Loss Equation 
(USLE): 

A = RKCPSL [22] 

where 

A = estimated soil loss (T/ac) 

R = rainfall factor 

K = soil credibility factor 

C = cropping management factor 

P = erosion control practice factor 

S = slope steepness factor 

L = slope length factor 

The S and L factors can each be expressed mathe- 
matically based on data analyses by Wischmeier and 
Smith (55). 
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n _ 



0.43 + 0.30s +0.043S 2 
6.613 



L = (X/72.6) 



05 _ 



8.521 



[23] 



[24] 



where 



s = percent slope 
X= slope length 

The USLE is widely used, and each of the factors has 
been determined for many soils and situations through- 
out the United States. It was originally derived to 
predict average annual soil losses. However, all factors 
can also be evaluated on a storm basis. Throughout this 
report, it is being used on a storm basis. In the USLE, 
the energy available for detachment is that derived only 
from raindrop impact. 

Williams (53) has modified the USLE to piedict 
sediment yield, G, at the watershed outlet. The rainfall 
energy term was replaced by an empirically derived term 
based only on flow characteristics for a particular storm 
event. 



G = 



KCPSL 



[25] 



where 



G = sediment yield (T/ac) 
Q ~ runoff volume (ac-ft) 
q p = peak flow rate (cfs) 
A r = drainage area (ac) 
a, = coefficients 

K, C, P, S, L = USLE factors 



The coefficients, a and 0, determined from 778 events 
on 18 watersheds in Texas and Nebraska, were 95 and 
0.56 respectively (53), The coefficients, calculated from 
data by Piest, et al. (37) from 103 events on two Iowa 
watersheds were found to be 23,6 and 0.65, respectively. 
The exponents are nearly equal, indicating that a change 
in the product, Qq_, has about the same effect. 
However, location differences are revealed more Signif- 
icantly by differences in a. 

The USLE and the Williams modification represent 
the two extremes with regard to energy required to 
initiate detachment, in that one utilizes rainfall energy 
entirely and the other, runoff entirely. Both rainfall and 
runoff have detachment potential and act simulta- 



neously at times during a rainstorm. The relative 
influence of each is dependent on both time and space in 
a watershed. 

At the beginning of a storm, before runoff begins, 
raindrop splash is detaching and transporting soil. This is 
the major energy source since runoff has not begun. 
When runoff begins, it also has potential that increases 
with the rising hydrograph. At the end of the storm, 
runoff is the detaching and transporting agent. 

With regard to position on a watershed, upslope 
erosion, where detachment is limiting, depends primarily 
on rainfall energy and slope gradient (54). Splash action 
and sheet flow transport the soil particles to flow 
concentrations or nils, Downslope, as runoff increases, 
its detachment and transpoit potential increases relative 
to rainfall. Finally, transport capacity becomes limiting 
and sediment yield becomes closely related to the 
hydraulic properties of the runoff and less related to 
rainfall characteristics (Wa). 

Foster(/0a) modified equation [1] to reflect both 
rainfall and runoff contributions. 



A = [aR + bcQq. 1 ' 3 ] KCPSL 



[26] 



where 



a,b = weighting parameters (a+b=l) 
c - equality coefficient 

R = rainfall factor 

Q- runoff vol. (in.) 
Q p = runoff rate (in./hr.) 
K, C, P, S, L - USLE factors 

The weighting parameters reflect the relative amounts of 
erosion caused by rainfall and runoff under unit condi- 
tions (L, S, C, P = 1). Limited experimental evidence 
(56, 28) indicates that the detachment is about evenly 
divided under these conditions (i.e., a-b=.5). The 
equality coefficient, c, can be estimated from rainulator 
plot data under unit conditions. 



c-(2A u /K-R)/Qq B 1 / 3 



[27] 



where 



A u = soil loss under unit condition. 



Analysis of 55 unit rainulator runs on 20 soils in 
Indiana 7 and Minnesota shows that c equals 29.3 or 



7 Foster, G. R., 1973. Personal communication. 
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approximately 30, Substituting these values of a, b, c 
into equation [2b] yields the basic equation for the 
erosion-deposition submodel. 



Foster and Wischmeier (10) have shown Uv.it the 
detachment rate at any point x is 



where 



A - EKCPSL 



E = 0.5R + ISQq. 



[28] 



1/3 



= 3A' 
x 9x 

or from equation [30] 
EKCPS (1.5 



185.58 



[31] 



[32] 



The coefficients, a, b, and c, may be dynamic in that 
they may change from time to time and place to place 
on a watershed during a storm. They may also be 
different during the detachment and transport phases of 
erosion. So might the other factors in the USLE. 
However, formulations of these relations- cannot be 
made at this time due to lack of data. Therefore, 
throughout the succeeding discussion, it is assumed that 
these coefficients and the USLE parameter remain 
constant at their indicated values through time and space 
and also for both phases of the erosion process. 



Detachment Component Development 

For a field of uniform slope, s, as shown in figure 10, 
equations [24] and [28] are combined to give: 

EKCPS X - 5 
A 8.521 

This can be extended to any point x on the slope so that 
EKCPS x ' 5 [29] 



A 



8.521 



When the soil loss is expressed as Ibs/foot of width, 
equation [29] becomes 



A' = 



EKCPS x 1 - 5 
185,58 



[30] 



Integration of equation [32] between two length .Hmlis 
yields the detached soil capacity for the length segment. 
So, for segment j of the complex slope shown in figure 
11 



[33] 



A dj = D x dx 



where A di = detached soil capacity for segment j or 






r 34 ] 



185.58 



where 



When this is extended to a compound slope having n 
segments, the cumulative detached soil capacity is 



A t - 



dM + A dj + A dj + 1 



dn 



This would equal the sediment yield per unit width if no 
deposition occurred along the slope length. Alterna- 
tively, the sediment yield equals the cumulative de- 
tached soil capacity provided the transport capacity Is 
not limiting, 



r 



W 




i-H 



+ 2 



Figure 10,-Field of uniform slope. 



Figure 11. -Complex slope profile, 



Transport Component Development 

For a uniform slope, let the USLE represent the 
transport capacity at any point x, so that 



E K CKPS 
185.58 



. . 
1 - 5 



[35] 



where 



A cx = transport capacity at x 



This assumes that C, P, and S remain constant and 
identical to those for detachment in the transport phase. 
The value of K should reflect the transportability of the 
material coming from upslope and the effect of soil 
properties on rilling and transport capacity on down- 
slope segments. Therefore, an effective value, K, is used 
which is the average K of the upslope segments. For a 
nonuniform slope, the transport capacity for any seg- 
ment is assumed to be that sediment load calculated 
using the total slope length for L and the slope steepness 
of the segment for calculation, S. Refeiring to figure 11, 
the transport capacity at x i+1 is 



CX>+1 






[36] 



185.58 

assuming a slope steepness of Sj throughout the lengtli 
Therefore, two cases can exist for slope segment j 

Case 1 



If A ( <A C , erosion occurs on slope segment j 
because the transport capacity exceeds the total de- 
tached soil. Therefore, the sediment yield 
G = A^i + +A d j_ 1 + A d j and the deposition M = 0. 



Case 2 



If 



A t >A c , deposition occurs on slope segment] 



because total detachment exceeds the transport ca- 
pacity, Therefore, the sediment yield G = A^^j and the 
deposition M - A t - A cx , . . 

Rilt-lnterrill Contributions 

The division of soil loss between rill and interrill areas 
is important for the chemical-sediment interface because 
of different methods of chemical application. For 



example, a chemical application tilled-in would have a 
different contribution from another applied as a spray. 

The conceptual basis for the nll-interrill contributions 
are separate laboratory studies conducted by Young (56) 
and Meyer, et al. {28} in which they observed the trends 
illustrated in figure 12. The top part of the figure shows 
that the interrill erosion rate builds up rapidly and then 
abruptly levels off, indicating that the rate remains 
constant downslope. The bottom part indicates that at 
some distance downslope rill erosion volume equals 
interrill erosion volume and, again, that interrill erosion 
increases linearly. Furthermore, when interrill and rill 
erosion volumes are equal, interrill erosion rate has 
already become constant. 

Foster 7 derived a relationship which approximates 
the ratio of 101 erosion to interrill erosion fora segment 
where detachment rate is at capacity. 



ft = 



j-x?.!) (0.043s?) 




[37] 



where s = percent slope and 

Kj/Kj = ratio of rill to interrill soil erodibility 

(rUlibility). 
Rough guidelines for the ratio, K r /Kj, are; 

K r /Kj = 2; "severe rilling" 

VK^l; "moderate rilling" 

K r /K[ = 0.5; "little evidence of rilling" 
In general, increases as slope length, steepness, rilli- 
bility and runoff increase. 

The quantity of sediment detached by rill erosion on 
a segment is 



[38] 



and the quantity detached by interrill erosion on a 

segment is 



[39] 



The equations for cumulative rill erosion and interrill 
erosion in a compartment are: 



m 
S 



m 



andIA= 



H 



[40] 



where m is the number of slope segments, 
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10 20 30 40 30 

DISTANCE DOWNSLOPE (ft) 

Figure 1 2.-Eroded soil apportioned between nil and iiiterrill 
sources for nil-susceptible Plot R-3 (From Meyer et. al. (28).) 

The interface to the chemical submodel requires that 
rill surface area and depth be specified. Limited experi- 
mental evidence obtained by Foster 7 on rainulator plots 
indicates that the percentage of the surface intercepted 
by nils at vanous cross-slope transects remains relatively 
constant from the top to' the bottom of the slope, 
although the number of nils becomes less in the down- 
slope direction. Based on these data and other data 
from construction site erosion, he postulates that the 
average rill width can be represented as follows: 



/Kr 

2.41- 



- 3 



[41] 



where 



w- average rill width (in./ft.) 

This indicates that the rill width depends on the 
rillibility of the soil, Kj/K r , and the slope steepness, s. As 
riUibilily and slope steepness increase, the average rill 
width increases. 

Average rill depth can be calculated by using the 
results of'equations {40] and [41] and utilizing the fact 
that 1 yard 3 of soil weighs about 1 ton assuming the 
bulk density of 1. 2. 



RD = 



324 (RA) 
wX 



[42] 



where RD = rill depth (in.) 
Eroded Material Particle Size Distribution 

It is important that not only the total amount of 
sediment be predicted, but also the composition of 
eroded mateiial should be known for the sediment- 
chemistry interface. This problem was approached utiliz- 
ing texture information from 56 Midwest soils. The 
original soil textures ranged from sand to silt loam to 
clay. The data consisted of measured soil loss and 
textural classification of both the soil matrix and the 
eroded material. This was done for each of two 
rainulator runs on four replications for each soil. 

For purposes of analysis, the following specific 
surfaces were assumed C/7): 

clay - 200m 2 /g 
silt - 40m 2 /g 
sand - 0.5 m 2 /g 

Therefore, the specific surface for a particular soil, i is 
SB; = 200 (%C1) + 40(%Si) + 0.5(%Sa) [43 J 



The clay ratio can also be calculated by 



CR ' = 



[44] 



Figures 13 and 14 show the graphs of the specific 
.urface and clay ratios of the soil matrix versus those for 
the eroded materials for the 56 Midwest soils. The 
regression equations are as follows: 



SS e 14.6 + 0.84 SS 



m 



and 



CR e = 0.021 + 1.08 CR 



-m 



[45] 



[46] 



where the subscripts m and e refer to soil matrix and 
eroded material, respectively. The coefficien ts ' of de- 
termination for equations [45] and [46] are 0.81 and 
0.82, respectively. Equations [43] -[45], together with 
the fact that the sum of the components equals 1, 
constitute the estimation procedure. This is a first 
approximation recognizing that many factors are not 
included. It is illustrated by the following example. 
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Example: Estimate' the mechanical composition of the 
eroded soil if the composition of the soil 
matrix is 

Cl m = 20% Si m = 40% Sa m = 40% 
From equation [43] 

SS m = 200(0.20) + 40(0.40) + 0.5 (0.40) =56 
and equation [44] 



= ' 25 
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Figure 13. -Plot of specific surface ; eroded vs. soil matrix. 



From equations [45] and [46] the estimated parameters 
for the eroded materials are: 

SS C = 14.6 -f- 0.84(56) = 62 
CR e = 0.021 + 1.08(0.25) = 0.29 

Therefoie, solving the following three equations for the 
eroded soil tcxtural components 

200 Cl c + 40Si c + 0.5 Sa e = 62 

Cle = 029 

Si e + Sa c 

Cl e + Si e + Sa c = 1 

show that the estimated eroded soil compositions are 

Cl e = 22% 
Si e = 43% 
Sa e = 35% 



and the enrichment coefficient is 
Cl e 70 



[47] 
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Figure 14. -Plot of clay ratios; eroded vs. soil matrix. 



Watershed Calculations 

The erosion-deposition submodel is applied to the 
watershed utilizing the transformation illustrated in 
figure 1. Each compartment may have a unique set of 
parameters required for the solution of the 
erosion-deposition equations. Calculations proceed from 
the top, downslope to the channel, continuing from one 
zone to the next. It is assumed that all soil delivered to 
the channel reaches the watershed outlet. 

The hydrologic submodel predicts the flow volumes 
and discharge rates from each zone but not from each 
compartment, Apportionment of flow volumes and peak 
discharge rates among the compartments within a zone 
were made using linear interpolation. The same tech- 
nique was used for adjacent compartments between 
zones. Because the zone is the smallest calculation untt 
permitted by the hydrologic submodel, all hydrologic 
input variables are identical for compartments traversing 
the slope, Input variable differences along the slope can 
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be handled by adding more zones to the hydrologic 
model. 

The Chemical Submodel for ACTMO 
"A Model of Chemical-Sail-Water Interactions" 

The objectives of the chemical submodel are to 
calculate concentrations of a chemical in runoff, the 
total amount of the chemical leaving the watershed, and 
the location of the remaining chemical, so that manage- 
ment practices can be selected to minimize undesirable 
movement. The submodel traces the movement of a 
single application of a chemical through and over the 
watershed. Cultivation is the only management practice 
explicitly accounted for since it redistributes the chemi- 
cal. All other practices such as contouring and terracing 
are implicitly handled through the hydrology and 
erosion submodels. Nitrate behavior was judged to be 
unique enough to warrant a separate submodel that uses 
many of the features in the general model, 

Simple chromatographic theory is used to describe 
the leaching of the chemical and to provide a method of 
calculating surface concentrations for runoff and ero- 
sion. Chemicals are degraded or lost, redistributed by 
cultivation, and nitrate is mineralized and taken up by 
the crop during the period of time between storms. Dur- 
ing the storm period the chemicals are transported by 
water and sediment. Adsorption characteristics of chem- 
icals range from highly adsorbed to no adsorption. Highly 
adsorbed chemicals, such as DDT and phosphate, have 
very low solution concentrations and are moved princi- 
pally by eroding the soil on which they are adsorbed. 
Other chemicals, such as nitrate, are not adsorbed and 
move with the water. Chemicals with intermediate 
degrees of adsorption are moved by both vehicles. 

The multitude of chemical reactions in the watershed 
are classified into three groups; adsorption or the 
distribution of the chemical between the solid and the 
solution phases, breakdown or loss of chemical, and the 
dispersion of the chemical as it moves through the soil. 
Each of these three groups of reactions is described by a 
simple relation which in many cases is only a crude 
approximation. More complete descriptions require in- 
formation not generally available. 

The compartment is the basic unit for computations. 

The hydrology submodel divides the watershed into 

ones white the erosion submodel creates routing tubes 

carry the cascading water and sediment across the 

**& to the channel. The compartment is the soil 

V of each zone within a routing tube shown in 

of the watershed transformation section. For 

3f chemical transport, the hydrology submodel 



is limited to a single ciop per zone and assumes the soils 
are similar within a zone. Although not essential, the 
layers of chemical interaction with the soil are assumed 
to be the same as the hydrologic layers. The English 
units, such as pound, acre, and inch, are used. A 
bibliography of chemical interaction and chemical trans- 
port (Appendix) has been compiled during the develop- 
ment of this submodel. These papers should be of 
assistance to anyone wishing to estimate the parameters 
for a given chemical and soil or those wishing to modify 
any concept used in this model. 

The Chemical Option 

Adsorption. The adsorption relation is assumed to be a 
linear isotherm, 



S = AC C 



[48] 



where S is the pounds of chemical adsorbed per pound 
of solid, C is the pounds of chemical per pound of 
solution, and AC is the adsorption coefficient. Batch 
equilibrations and column leaching studies t (Appendix) 
provide information on the magnitude of the adsorption 
coefficient for many soils. The model can be easily 
changed to include nonlinear relations, but this requires 
another parameter. 

Rates of chemical adsorption and desorption (26} 
were considered to introduce more complexity than is 
merited for this simple model. Nonequilibrium condi- 
tions cause the chemical distributions to be asymmet- 
rical .and are probably more prevalent in the real world 
than not. However, it would appear that this effect is 
usually not any larger than the normal variability 
encountered in the field (40, 41}. 

Breakdown. There are numerous papers (Appendix) on 
the breakdown, degradation, loss by volatility, and 
persistence of chemicals in both laboratory and field 
conditions. A simple first order rate equation was chosen 
to represent these many processes. The chemical remain- 
ing, A #/a, is related to the initial amount, AA #/a, by 
an exponential function (EXP) of time, TI days, 



A = AA EXP (-BC - TI) 



[49] 



The rate coefficient, BC/day, is dependent upon the 
environment around the chemical (.5). There isn't 
enough quantitative information available at this time to 
make the necessary corrections, but the model has been 
designed to accommodate the calculations when the 
relations are developed. The average temperature and 
soil moisture are available from the hydrology submodel. 



The type of application, banding or broadcasting, also 
changes the degradation rate (<5) as can repeated applica- 
tions of the same chemical (16), Such factors have been 
neglected in this first approximation. -Large amounts of 
pesticides can be lost by volatilization during application 
and, theiefore, the amount at the start of the simulation 
is not the amount leaving the applicator, but rather the 
amount that actually reaches the soil. Volatilization loss 
after the chemical reaches the soil is usually significant 
only if othei losses are small. Volatilization is a complex 
process involving a number of factors, and this first 
version assumes that these losses will be reflected in the 
apparent breakdown coefficient. A detailed discussion of 
pesticide volatilization is given by Spencer et al. (43). 
Dispersion. When a chemical is carried through the soil 
by flowing water, there is a broadening of the chemical 
band and a decrease in the peak concentration. Two 
processes are in operation: diffusion and dispersion. 
Diffusion is the movement of the chemical from areas of 
high concentration to areas of low concentration, 
Hydrodynamic dispersion is the result of unequal water 
flow through the different sized pores and cracks in the 
soil and should reflect the aggregate size (36) and the 
ratio of porosity to permeability (4). The dispersion 
coefficient has the same units as the diffusion coefficient 
(area/time) and is proportional to the waterflow 
velocity. For the flow rates during infiltration and 
drainage, the dispersion coefficient is considerably larger 
than the diffusion: coefficient. A more useful factor for 
the equations to be used is the ratio of the dispersion 
coefficients and the average flow velocity. For the sake 
of clarity this is called the dispersion distribution 
coefficient, DC, with units of length. 

Equations involving dispersion are mathematically 
similar to diffusion equations where many solutions are 
available (24). Some experimentally measured values of 
DC from column studies range from .03 inch with 
sieved soil (7) to about 5 inches with undisturbed soil 
columns (14). Little field information is available for 
tins parameter. 

Movement. There are basically two ways to calculate 
the movement of a chemical through the soil. The first is 
the finite element method in- which the soil is divided 
into thin layers of soil (9, 1 2). The second is to assume a 
simple continuous chromatographic process. In the first 
method the chemical is moved from one layer to another 
and complete mixing with uniform distribution within a 
layer is assumed, Two problems arise with this approach. 
When the chemical is applied in a narrow band, a large 
number of layers are needed to accurately trace the 
chemical movement. The second problem is to define 
the thickness of the soil surf ace -run off interface. The 



ratio of the watershed surface area to the specific surface 
area (area/area per unit weight of soil) is one estimate of 
the amount of soil in this interface {23}, but this appears 
to be a considerable underestimate. 

An effort was made to circumvent these problems by 
assuming an idealized bell-shaped distribution or plug 
distribution for use in a simple chromatographic process. 
This assumption of a known mathematical equation 
provides a means to calculate the concent lation at any 
depth including the surface of the soil. It appears that 
such a proceduie is adequate when there Is no adsorp- 
tion or the adsorption isotherm is linear. This holds for 
nitrate, phosphate, and many pesticides, but not for 
cation exchange involved in salinity movement (39). 

Water moves the chemical through the soil along 
several paths. The infiltrating water tends to move the 
chemical vertically and sometimes laterally through the 
soil. Runoff water carries some chemical with it and also 
causes erosion which can carry additional amounts of 
chemicals. Table 6 gives an outline of the calculations 
made in linking chemical characteristics with the move- 
ment of the carrier. Two types of chemical applications 
are provided for on the surface or mixed into the soil. 
The amount applied is assumed to be corrected for 
application losses. The surface application includes the 
amount on the plant surface, if any, that is readily 
washed off. A major gap exists in the understanding of 
factors controlling the rate of washoff from the plant 
surfaces. Mixing the chemical into the surface soil is 
assumed to result in an ideal uniform plug-type distribu- 
tion even though this may not always be true. 

Table 6,-Outline of steps for tracing chemical 
movement 

1. Amount of loss before the storm, 

2. Movement by infiltration. 

3. Concentration at the surface. 

4. Runoff. 

5. Intcrrill erosion. 

6. Rill erosion. 

7. Accumulation and transport of water, sediment, and 
cliemical along routing tubes. 

8. Deposition, 

9. Lateral outflow. 

10, Mixing of water, sediment, and chemical from otlier 
tubes. 

1 1. Location of the remaining chemical 

Losses before storms, step 1. The chemical is lost 
according to equation [49] for the period of time from 
application until the first storm and between storms 
thereafter. The time-average soil moisture and tempera- 
ture are available for this period and could be used to 
modify the amount of breakdown or loss. An extra 
coefficient for the photodecomposition rates at the 
surface is used m the array of breakdown coefficients. 
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Infiltration, step 2.A great deal of information has 
been developed about the infiltration process (7 J), Since 
it was decided to make this first attempt a very simple 
model, the leaching process was divided into two stages. 
The first stage is prior to runoff when nonadsorbed 
chemicals like nitrate and halides move with the water 
front into dry soil (35}, This distance is approximately 
the amount of infiltration in inches, IN, divided by the 
volumetric field capacity, FC. 



D = IN/CFC+BD-AC) 



[50] 



The term (BD-AC) accounts for the chemicals lagging 
behind the waterfront when the chemical is distributed 
between the water and the soil of bulk density BD (13). 
Some anions like nitrate are excluded from some of the 
fine pores in the soil, and BD'AC can be given a negative 
value corresponding to the water volume from which the 
antons are excluded {42). The term field capacity is used 
as an operational definition of the amount of water 
remaining after saturated flow has ceased. It is recog- 
nized that tliis is physically imprecise and difficult to 
evaluate. 

In this model the second stage of leaching is for the 
infiltration that occurs during runoff. There is consider- 
able data showing that continuous percolation is less 
efficient in leaching salts than intermittent flow. This 
inefficiency is probably a result of relatively more water 
moving through large pores than small ones. An ex- 
amination of the data in Miller et al. (29) shows that the 
inefficiency is proportional to the water content above 
field capacity. Therefore, the leaching efficiency used is 
the ratio of field capacity to the average moisture 
content during this period. The average moisture content 
is calculated as halfway between field capacity and 
saturation since the process of reaching saturation and 
the return to field capacity accounts for the bulk of the 
water movement. The depth of chemical movement, D 
inches, is calculated by 



[51] 



FC AM+BD-AC 



where IR is the inches of infiltration during runoff, AM 
is the average moisture content, and the other variables 
have been defined. 

Experimental evidence has been obtained to verify 
some of these relations. Figure IS illustrates the be- 
havior of surface-applied salts during infiltration on soils 
>f different moisture contents (unpublished data from 
)urant, Oklahoma), Lithium bromide salt was applied to 
IB surface of columns of a Fort Collins silt loam soil 
then sprinkled with one inch of water at the rate of 



1/4 in/hr. When the soil was air dry, the bromide, which 
is leached similai to nitrate, moved with the water front 
and was found at a depth of 4 inches in a sharp band. 
When the soil was wet at nearly field capacity, the 
simulated rain soon ponded and suction was used to 
draw the water through the soil. In this case there is 3 
loss of leaching efficiency and a more dispersed band is 
found at a depth of 2,5 in. Lithium is a slightly adsorbed 
chemical with an adsorption coefficient of about three 
in these experimental conditions and thus is moved only 
about one-tenth of the bromide distance. This distance is 
too small to show any effect of the wet or dry 
treatment. 

The concept being used for leaching efficiency and 
dispersion was tested for field conditions with published 
data of Miller et al. (29), Potassium chloride was applied 
to small plots and leached using three irrigation treat- 
ments: continuous ponding, successive 6-inch incre- 
ments, and successive 2-inch increments. The chloride 
content of the water at different depths was monitored 
with porous cups and soil water by tensiometers. The 
dispersion distribution coefficient was calculated from 
all the peak concentrations relative to the distance they 
had moved. The calculated peak location, height, and 
width are compared with five sets of observed data in 
figure 16. The fit with othei observations is equally good 
which suggests that the relations used should be ade- 
quate for the first approximation model being de- 
veloped. 
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Figurc 15. -The movement of lithium bromide through a column 
of Ft. ColUns silt loam by one inch of simulated rain when 
the antecedent moisture was air dry or at field capacity. 
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Concentration at different depths and times, step 
3.-The combination of simple chroniatographic be- 
havior with dispersion and a Gaussian or bell-shaped 
distribution gives rise to equation [52] from which the 
solution concentration, C ppm, at any depth (X inches) 
can be calculated. 



A-U 



V4TT-DC-D 



EXP - 




[52] 



A is the amount of chemical and U is a units factor 
equal to 4.42/(SM+BD*AC) when A is in pounds per 
acre and C is in ppm at a volumetric moisture content of 
SM. When X = 0, the equation gives the concentration 
at the surface and for X = D the peak concentration in 
the soil is calculated. 

When the chemical is mixed into the soil AP inches 
deep and there Is infiltration that moves the plug 
downward D inches, the boundaries of the plug are 
dispersed in a sigmoid pattern (33). The concentration at 
X inches from the middle of the plug is given by (24): 



C(X) - 1/2 CI- 



ERF 



ERF 




\/4'DC 



1/2AP-X_ 
\/4'DC*D 



[53] 



where ERF is the error function and CI is the initial 
uniform concentration in ppm: 



CI = A/(.3333 AC + .2266 SM)AP 



[54] 



with .3333 million pounds of soil per acre inch and 
.2266 million pounds of water per acre inch. Equation 
[53] is used to calculate the concentration at the 
surface, X D + 1/2 AP, and to check if the plug 
distribution can be treated as a bell-shaped distribution. 
When the concentration at the middle of the plug X - 
has been depleted more than five percent by the 
dispersion effects, then the distribution is assumed to be 
bell shaped. 

After a number of storms the concentrations In the 
solution and on the solid at the surface may be reduced 
below a predetermined critical level, and the runoff and 
erosion calculations can be bypassed. 

Runoff, step 4.-Affer the initial infiltration when the 
water begins to pond on the surface and runoff 
commences, it is assumed that the rainfall -equilibrates 
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Figure 16. -The calculated and observed solution concentrations 
of chloride with depth under two irrigation treatments 
(Panache clay loam, California). 

with the solid phase on the soil surface. Then the 
concentration in the water as it moves away from the 
surface is the same whether it moves vertically as 
infiltration or laterally as runoff. Under these conditions 
a projected depth of leacliing can be calculated with 
equation [51] using the sum of infiltration and runoff. 
The concentration 'in the solution al the surface, and 
thus in the runoff, can be calculated with the projected 
depth in equations [52] or [S3] < 

The loss of the chemical to runoff should be 
calculated by integrating the surface concentration as 
the depth increases to the projected depth. No simple 
way has been found to perform this integration so an 
average concentration is calculated from equations [52] 
or [53] using an average distance between the depth at 
the start of runoff and the final projected depth. 
Examination of the behavior of these equations Indicates 
that there will be an underestimation of the mean 
concentration when the runoff and infiltration remain 
constant. However, the infiltration is usually higher at 
the beginning of the storm, carrying a larger fraclion of 
the high concentration into the soil, which reduces the 
underestimation. The total amount of chemical lost in 
runoff is calculated as the product of the mean 
concentration during the runoff and the amount of 
runoff. 
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Interrill erosion, step 5 .-Sediment from erosion has 
higher concentrations of clay, organic matter, nitrogen, 
and phosphorus than the original soil. The enrichment 
coefficient is the ratio of the concentration in the 
sediment to the concentration in the original soil. The 
highest enrichment coefficients tend to occur with the 
smallest amounts of erosion. Increasing amounts of 
erosion tend to bring the ratio closer to unity [55]. 
This suggests that the sediment from nil erosion has the 
same composition as the soil and that the enrichment is 
(he contribution of interrill erosion which is pre- 
dominantly silt and clay. Enrollment coefficients of two 
are quite common for clay, organic matter, etc., from 
erosion plots. 

The erosion mode) supplies estimates of the clay 
enrichment coefficients, EC, the pounds of interrill and 
rill erosion, GI and GR, for each compartment in the 
watershed, and the clay fraction, CF, and silt fraction, 
SF, of the surface soil in each zone. The pounds of clay 
in the interrill erosion, PC, is the total clay in the 
sediment minus the clay from rill erosion: 



PC = (CF EC) (GI + GR)-CF'GR 



[55] 



and the pounds of silt, PS, is the remainder. 

There are numerous references, see Bailey and 
White (.A), to the fact that the different size fractions of 
the soil adsorb chemicals to different degrees. The 
amount of chemical per unit of soil at the surface conies 
from the mean concentration during runoff and equa- 
tion [48] . The fraction of this that is adsorbed to the 
clay is PF CF/(PF'CF+SF) where PF is the preference 
factor of the chemical for the clay relative to the silt 
when there is none adsorbed on the sand. Therefore, the 
pounds of chemical per pound of clay, SC, is: 



SC = 1- PF ' CF 

CF PF CF + SF 



[56] 



and the remainder of the chemical is associated with the 
silt. Similar calculations with an additional preference 
factor can be used when the chemical is associated with 
the sand fraction. 

The tolal amount of chemical removed in interrill 
erosion is the sum of the product of concentration and 
amount for the clay and silt. 



AI - SC-PC + SS-PS 



[57] 



Rill erosion, step 6.-The basic logic used for the loss of 

-jhemicats in rill erosion is that all the chemical within 

he rUl volume is lost. The area of the rills, EA, as a 

action of the treated area and the average depth of the 



rills, ED, is provided by the eiosion model. Rills , 
probably wider at the soil surface, but for simplicit; 
rectangular rill shape is assumed. The method 
estimating the width and depth of the rills is cntde 
this time and doesn't merit a more complex sha 
consideration. The amount of chemical within the 
volume is the product of the amount of chemical in t 
compartment, the area of the rills, and the fraction 
the vertical distribution intercepted by the rill. TJ 
fraction, Yl, is the integral of a unit distribution frc 
the surface to the bottom of the rill' 

ED 

Yl = f f(x)dx [51 



When the chemical has been applied to the surfai 
and leached D inches into the soil, then f (X) is a be 
shaped curve about D and the integral is 



Yl = 1/2- 



ERF 




ERF 



ED-D 



[55 



where the minus sign is used when the argument 
negative. A band width of two standard deviation 

2* V4'DOb, is used to determine how much of tli 
chemical distribution the rill depth intercepts. When th 
chemical is applied by mixing AP inches into the soil an 
then leached D inches, a plug with sharp boundaries 
f(x) = constant, between D and D + AP is assumed. Thl 
creates three limits: 



a. When ED is less than D then 
Yl = 



[60 



b. When ED is greater than D but less than D + Al 
then 

Yl = ED/(D + AP) [61] 

c. When ED is greater than D + AP then 

Yl - 1 [62] 

Cultivation that reaches the chemical is assumed to 
create a plug distribution and destroy the previous rills, 
Without cultivation the rill erosion is assumed to occur 
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in the old rills, and the effect of the combined nil depth 
must be used in setting the integral Yl. Summarizing, 
(he amount of chemical lost in rill erosion is 



AR = A'EA-Yl 



[63] 



Transport to the channel, step 7. Runoff water and 
sediment from one zone can cascade over other zones. 
The transformation of watershed into a series of tubes 
was used to describe this feature. The amounts of runoff 
water, sediments, and any chemical associated with them 
arc accumulated for all the compartments in a tube. 
Currently no allowance is made for loss of the chemical 
over untreated ground by infiltration or adsorption on 
the soil surface. One estimate of the amount lost by 
infiltration is the product of the infiltration during 
runoff and the concentration equilibrated with the water 
and sediment of runoff in that compartment. In most 
cases these effects are considered to be small. Since the 
soils could be different in their adsorption and pref- 
erence characteristics as well as texture, the amounts of 
clay and silt arc calculated. The absorption and pref- 
erence factors are weighted for the amount of sediment 
from each zone, and the chemical distribution between 
the sediment and water is calculated for the outflow of 
each compartment. 

Deposition, step 8. In some compartments the flow of 
runoff water will be insufficient to maintain all the 
sediment in suspension -and the excess sediment will be 
deposited. The hydrodynamic characteristics of the 
sediment cause the larger, heavier particles to settle out 
first. In most cases these heavier particles are sand grains 
containing a negligible amount of chemical. It is recog- 
nized that sometimes water-stable aggregates and 
organic-coated grains may cause a significant fraction of 
the chemical load to be deposited with the sand fraction. 
The chemical is distributed according to the weighted 
adsorption coefficient between water and sediment 
before the deposition. The amount of deposition is 
satisfied by depositing the sand, silt, and clay in that 
order. The weighted preference factor is used to cal- 
culate the deposition of the chemical with the different 
size fractions. 

Lateral outflow, step 9. The hydrology model provides 
estimates of the vertical and lateral flow rates for each 
layer in each zone. The compartments at the end of 
tubes are next to the channel, and chemicals in these 
compartments can be carried into the channel by lateral 
flow. Water is also contributed to the outflow of each 
Lube. 

The amount of lateral waterfiow Is the ratio of the 
lateral flow rate to the total flow rate times the amount 



of drainage The amount of drainage from each layer is 
the amount of infiltration to each layer minus the 
amount required to reach field capacity, Since the 
chemical could be distributed between two layers, the 
fraction in each layer must be calculated. For a bell- 
shaped distribution the fraction above the boundary is 
given by equation [58] where ED is replaced by the 
distance to the boundary. The fraction for a plug 
distribution is the difference between the boundary and 
the top of the plug divided by tiie thickness of the plug. 
The chemical outflow from each layer is then the 
amount of chemical in solution times the outflow 
divided by the amount of water in the soil. 

Mixing for watershed outflow, step lO.-The outflow 
from a watershed is not usually a "well-mixed system. 
Pairs of samples taken at the same time can differ in 
concentration from 9 to 33 percent (7 7). Some of this 
variability can be attributed to the cascading effect 
symbolized by the tube transformation in this model. 
Thus one estimate of the peak solution concentration in 
the outflow is from the average solution concentrations 
in each tube outflow. The combination of all the tubes 
gives the watershed outflow with a mean solution 
concentration calculated by distributing all the chemical 
between all the sediment and all the water. 

Location of the remaining chemical, step 11. The 
presence of a lateral component'of waterflow results in 
a diagonal movement of the chemical through the soil 
profile. The total distance moved has been calculated 
from the infiltration. The vertical and horizontal com- 
ponents are calculated from each increment of move- 
ment using the ratio of lateral to vertical flow rates 
weighted for the fraction of the chemical in each layer. 
The peak solution concentration is calculated from 
equations [52] or [54], The deposited chemical is 
added to the chemical already in that compartment or 
creates a new treated area if there is no previous 
chemical. A plug distribution is assumed with a thickness 
of 1 inch plus the depth of the chemical layer already 
there. 

The logic for cultivation effects.-If the. depth of 
cultivation is deeper than half the chemical, it is assumed 
that a plug distribution is created from the surface to the 
depth of the cultivation, 

Iteration.-The preceding calculations are performed for 
each compartment in the watershed although some steps 
can be omitted if no chemical is present. After the 
calculations are completed for the entire watershed, the 
results are totaled and the status of the system is printed 
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out The program then returns to the beginning foi 
another storm event. 

The Nitrate Option 

Nitrate is derived in the soil not only by fertilizer 
application, but also by the mineralization process of 
micro-organisms. Nitrogen in organic matter is converted 
to ammonium, then to nitrite, and finally nitrate. Nitrite 
seldom occurs in significant concentiations and while 
ammonium occurs at significant concentrations more 
frequently, it is moderately adsorbed to the soil. Thus 
the nitrogen problem is usually restricted to nitrate 
which is nonadsorbed and leached easily. Plant uptake, 
denitnfication.and immobilization are other mechanisms 
for nitrate removal. 

Since erosion can remove adsorbed organic forms of 
nitrogen and ammonium, but does not significantly 
affect the nitrate movement, the erosion model was 
bypassed and the storm file, from the hydrology model 
normally used by the erosion model, is used directly. 
The elimination of the erosion model also eliminates the 
need to consider compartments within the hydrologic 
zones. Normal mineralization means that the nitrate 
movement in all zones has to be calculated, whether the 
zones received fertilizer or not. The starting date for the 
simulation is the date of the first application of 
fertilizer. The model provides the option for a split 
application of fertilizer to be applied on a second date. 
Table 7 presents an outline of steps used in calculating 
nitrate movement. 

Mineralization, step 1. -Organic nitrogen is mineralized 
to nitrate according to a first order rate equation (45), 
The pounds per acre of nitrate, AN, from mineralization 
during a period of TI days is; 



AN = PN[1.~EXP(RK-TI)] 



[64] 



where PN is the pounds/acre of potentially mineralizable 
nitrogen in the top layer of soil, and RK in reciprocal 
days is the mineralization rate coefficient. The rate 
coefficient is sensitive to temperatures (46) and mois- 
ture (47). The Arrhenius equation relating the rate 
coefficient to TM, the average absolute temperature 
(273 + degrees Celcius) during the period, is 



ln(RK) = 15.807 - 6350./TM 



[65] 



The optimum moisture for mineralization appears to be 
near field capacity, and the rate coefficient' decreases 
linearly as the moisture content decreases: 



WC = SW/FC 



[66] 



Table 7.-Outline oj steps for tracing nitrate movement 

1. Mineralization, temperature, and moisture corrections. 

2. Plant uptake. 

3. Fertilizer and cultivation, 

4 Vertical movement and lateral outflow. 
5. Iteration and band completion. 



where WC is the correction for the soil moisture content, 
SW, averaged over the period of mineralization. When 
the moisture content is above field capacity, the conver- 
sion of ammonium to nitrate is inhibited but the 
conversion of organic nitrogen to ammonium is not. 
Short peiiods of ponding are incorporated in the model 
through the second infiltration stage. Any ammonium 
buildup is dissipated quickly to nitrate when aerobic 
condjtions return. 

Ponding does accentuate denitnfication which is not 
accounted for in this model. Considerable study was 
given to this process before it was decided to neglect the 
piocess In this version. The presence of a carbon source 
is necessary for the conversion of nitrate to gases, and 
this would require additional modelling of soluble 
carbon compounds, The same problem was faced and 
the same decision was reached for the process of 
immobilization, the conversion of nitrate into organic 
nitrogen compounds. These are areas for future research. 

Plant uptake, step 2.-The uptake of nitrate.by the crop 
is the mechanism of greatest nitrate loss during the time 
when plants are actively growing. Unpublished work 
indicates that when the accumulated mineralization plus 
fertilizer does not greatly exceed the amount of nitrogen 
for optimum plant growth (44) only a small amount of 
nitrate remains in the root zone at plant maturity. 

Laboratory studies (51) show that nitrogen uptake 
increases with the increasing concentration of nitrate in 
solution. Mass flow is the predominant mechanism for 
moving nitrate through the soil to the plant root (2). 
Thus nitrogen uptake should be closely connected to 
evapotranspiration, Evapotranspiration is also related to 
the amount of adsorbing root surface and reflects the 
growth rate of the plant. There is little growth, ET or 
nitrate uptake under cold or drought conditions (52). 
Based on these considerations the nitrate taken up by 
the plant, UP pounds/acre, during the period between 
storms is: 



UP * AU ET/(SW WD) 



where AU is the pounds of nitrate per acre available for 
uptake, the initial content plus half the mineralized 
amount; ET Is the inches of evapotranspiration from 
each layer weighted for the distribution of the nitrate 
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within layers imd the piocluct of SW and WD gives the 
inches ofwatei in the soil 

Fertilization and cultivation, step 3. -Fertilizer added 
before a storm is not considered available for uptake 
because uptake horn split applications of nitrogen is 
small if theie isn't adequate moisture, The applied 
fertilizer is available for leaching dining the storm and 
for uptake after the storm. Cultivation is not considered 
directly in the nitrate model, but its effects on water 
movement are accounted for in the hydrology model. 

Movement, step 4. -The same logic is used for the 
leaching of .nitrate as is used in the general chemical 
nibdcl. Infiltration causes vertical movement and lateral 
outflow which is added to the surface runoff. There is 
no nitrate cairied in the suifacc runoff except as 
contributed by return How, and no loss by erosion or 
deposition is calculated. The model could be modified to 
handle a surface application of nitiogcn fertilizer using a 



dissolution rate rather than assuming complete dissolu- 
tion in the initial infiltration 

Iteration and band completion, step 5 .-The preceding 
calculations are repeated for each zone to provide the 
watershed output. All the calculations are repeated for 
each storm event causing the band of nitrate to move 
downward. Since mineralization is assumed to occur in 
the top layer, the mineralization subsequent to each 
storm in effect creates a layer that is thicker by the 
amount of movement. When the mineralization rate is 
slow enough, its contribution to the band of nitrate can 
be bypassed and this defines the location of the top of 
the band. This decision is made when the average 
temperature is 5 C or less. The temperature relation 
indicates that the mineralization rate at 5 C is about 10 
percent of optimum. 

Any infiltration during this cool weather {fall, winter 
and spring), causes the band to move further. No explicit 
provisions have been made for frozen conditions. These 
arc presumably handled by the hydrology model. 



TESTS FOR VALIDATION 



Carbofuran Experiment at Coshocton, Ohio 

Complete data sets of chemical, erosional, and hydro- 
logical information are very scarce, The only complete 
set we have been able to compile is for a 2-year 
experiment with the pesticide carbofuran on watershed 
No, 113 at the North Appalachian Experimental Water- 
shed, Coshocton, Ohio. The data set was developed from 
the manuscript reporting the lesults(tf) and supple- 
mented with information from the authors and the staff 
nt the watershed. 



Input to tho Hydrological Submodel 

Tables 8 and 9 give the input format for the 
parameters and data. The l'.5-acre watershed was divided 
into three equal zones as shown in figure 17. The upper 
zone had an average length of 80 ft and slope of 3 
percent. The middle zone was 110 ft with a 10-percent 
slope and the lower zone was 190 ft with 20-percent 
slope. While there is subsurface movement of water In 
Uiis watershed, none of it comes to the surface above the 
measuring weir. Corn was grown on the whole watershed 
the 2 years of the experiment, 1971 and 1972. To 
initialize the hydrologic conditions, the hydrology sub- 
model was operated for the 2 years preceding the 
experiment when the watershed was in meadow. Tem- 
perature and pan evaporation were entered as weekly 



averages of the daily values, mean degrees Fahrenheit 
and inches. 

The rainfall input was breakpoint values. Each ob- 
servation contained the date, military time, and the 
inches of precipitation since the previous reading. The 
time increment between observations did not exceed 2 
months when there was no rainfall and 30 minutes 
during storm. 

Input to the Erosion Submodel 

The information about each compartment needed by 
the erosion submodel is given in table 10. 

In watershed 113 all the overland flow cascades to 
the next zone so there is only one stream tube or path. 
The hydrologic variables required are the runoff volume, 
peak runoff rate, and the rainfall intensity for each zone 
in each storm, The slope length and steepness are the 
same as used in the hydrology. In addition, the soi] 
texture, 15 percent clay and 60 percent silt, the 
erodibility factor, 0.48, the crop factor for corn chang- 
ing from 0.29 to 0.43 with season, and the managemen' 
factor of 0.5, 0.6, and 0.9 for these zones are used, 

Input to the Chemical Submodel 

Table 11 lists the information required by tl 
chemical submodel for each compartment. Only OIK 
layer per compartment is necessary because of the 
limited leaching, 
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190 FT. 
ZONE HE 20% 




0.64 m. fell during this period and was treated as initial 
infiltration. Dunng the 25 days between the second and 
third storms, 2.49 in. of ram fell and was treated as part 
of the initial infiltration. Eight other storms with less 
than 0.08 in. runoff were measured, but not predicted by 
the hydrology submodel before the November?, 1972, 
stoim. 

Table IO,~Erosion submodel input data for each 

compartment 

Constant watershed characteristics 

1. Soil credibility factor 

2. Slope steepness, percent 

3. Length of overland flow, ft 

4. Width of compartment, ft 

5. Texture 

Variable watershed characteristics 

6. Crop factor 

7. Management factor 



Figure 17. -"Watershed 113 divided into three hydro logical zones. Storm characteristics 



From the hydrology submodel the values of the 
initial soil moisture, infiltration, runoff > and subsurface 
movement ate needed. The erosion submodel supplies 
values of interrill and riil erosion, area and depth of nils, 
deposition, and composition of the sediment. The 
chemical characteristics used are given in table 12. In 
1971 the carbofuran was broadcast on the soil surface 
and disced 3 inches into the soil before the corn was 
planted. In 1972 the carbofuran was banded 2 inches 
deep with the corn seed. The adsorption coefficient was 
calculated from the concentration on the sediment and 
in the solution while the breakdown coefficient was 
reported in the manuscript. The banding method in- 
creased the adsorption coefficient threefold and de- 
creased the breakdown coefficient twofold. The disper- 
sion distribution coefficient and the preference factor 
were chosen arbitrarily as "reasonable" values after 
consulting the literature. 

Results 

The predicted values of five variables are compared 
with measured values for five storms in table 13. In 1971 
rain started in the evening following application and 
continued the next morning to give the only major 
storm that season. Five-other storms with less than 0.05 
in. of runoff were measured but not predicted by the 
hydrology model. 

In 1972, 26 days passed from the time of application 
until the first storm runoff was measured. Rain totaling 



8. Rainfall factor from storm intensity 

9. Volume of runoff, in 

10. Peak rate of runoff, in/ltr 



Table 11. -Information required for the chemical 
submodel 

A. Chemical 

1. Date of application 

2. Application rate, Ibs/a 

3. Depth of mixing, in 

4. Adsorption coefficient, ratio of concentrations 
5- Breakdown coefficient, day" 

6. Relative preference for size fractions 

B. Physical 

7. Dispersion distribution factor, in 

8. Field capacity, % 

9. Poiosity, % 

10. Texture of soil 

C. Hydrological 

11. Date of storm 

12. Infiltration before runoff, in 

13. Infiltration during runoff, m 

14. Runoff, in 

15. Ratio of lateral to vertical flow 

16. Soil moisture at start of storm 

D. Erosional 

17. Interrill erosion, tons 

18. Rill erosion, tons 

19. Fraction of area in rills 

20. Depth of rills, in 

21. Deposition, tons 

22. Texture of sediment 
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Table 12 -Chaiactenstics of cat bofuran experiments 
on Coshocton Watershed 1 13 



Date of application 



May 5, 1971 May 18, 1972 



Amount applied, Ib/a 


4.83 


2.77 


Application method 


Broadcast 


Handed 


Depth of mixing, in 


3 


2 


Adsorption coefficient 


1 4 


4.1 


Bicdkdown coeflicicnt, day 


0015 


0.0074 


Dispersion distribution coefficient, 






m 


.5 


5 


Prefereiue factor 


10 


10 



The predicted chemical loss and concentration arc 
high. Only preliminaiy efforts have been made so far to 
evaluate the cause of the disciepancies fiom the three 
models. Additional results from other chemicals, storms, 
and wateisheds are needed before a thorough evaluation 
can be made 

A moie detailed analysis of the 1971 storm is given in 
table 14. About 64 percent of the chemical loss is 
removed by runoff water in this storm. A similar analysis 
of the 72 storms shows a decrease in the fraction of the 
chemical loss fiom runoff water as the chemical is 



Table 1 3. -Measured and predicted water, sediment, and carbofuran losses from 
Coshocton Watershed 113 for 5 storms 



Item 



Storm date 



5-6-71 



6-13-72 6-15-72 7-10-72 11-7-72 



Runoff, in 












Metis 


0.21 


0.24 


0.41 


0.20 


0.22 


Prcd 


.21 


.46 


.27 


.001 


.15 


Peak flow, cfs 












Meas 


.28 


.85 


3.54 


4.24 


... 


Prcd 


.23 


1.15 


.83 


.01 


27 


Erosion, tons/a 












Meas, 


1.73 


2.12 


4.49 


.69 


,38 


Pred 


4.27 


9.53 


4.80 


.80 


3.60 


Carbofuran, Ibs 












Meas 


,036 


.021 


.046 


.004 


.001 


Prcd 


.071 


.15 


.104 


.04 


.013 


Concentration, ppm 












Meas 


.47 


.19 


.22 


.05 


.02 


Pred 


.78 


.56 


.70 


4.24 


.14 



Table 14,-The movement of water, sediment, and carbofuran in each compartment for [he 1 91 '1 storm 



Item 


Compartment 


Total 


1 2 


watershed 
3 



Initial infiltration, in . . . 
Remaining infiltration, in 

Runoff water, in 

Runoff chemical, Ib 

Intenill soil, tons . , 

Intemll chemical, Ib. . , . , 

Rill soil, tons 

Rill chemical, Ib 

Leaching depth, in 

Remaining chemical, Ib/a . 



1.64 

.99 

,38 

.0265 

.1 

.0002 



.85 
4.7 



1.55 

.84 

.62 

.0454 

,1 

.0002 
1.1 

.0046 

,87 
4.6 



1.99 
1.38 

-.36 

-.0268 

.5 

.0003 
4.9 

.02 
1.03 
4.8 



.0451 
.0007 
.0246 
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leached deeper into the soil The negative values for 
runoff water and chemical in compartment three oc- 
curred because the infiltration in this compartment 
exceeded the rainfall and part of the cascading water 
from the second compartment was used. 

Erosion Studies at Treynor, Iowa 

Since additional complete chemical, erosional, and 
hydrological data sets are not presently available, otjier 
data sets have been used to evaluate parts of the model. 
The hydrology submodel has been evaluated with a 
number of data sets (20) and will not be discussed here. 
A set of hydrological and erosional data from watershed 
2 at Treynor, Iowa, was obtained from the staff of the 
North Central Watershed Research Center, Columbia, 
Mo., and used to evaluate the combination of the hy- 
drological and erosional submodels. 



Input to the Hydrological Submodel 

A map of the 82.8-acre watershed is shown in figure 
18. The watershed was divided into three zones. 200ft; 
166 ft, and 133 ft in average overland flow length. It was 
estimated that twenty percent of the upper zone 
cascaded to the middle zone while 15 percent of that 
zone cascaded onto the lower zone. In 1965 corn was 
grown, and in the 10-day period from June 25 to July 4 
six major runoff events occurred. To initialize the 
hydrological conditions the hydrology submodel was 
started the previous year when corn was also grown. 



Input to the Erosional Submodel 

To compute the erosion variables the watershed must 
be transformed. The transformation creates five stream 
tubes or paths and seven compartments as shown in 
figure 19. For this silt loam soil the texture is 15 percent 
clay and 65 percent silt, the credibility factor is 0.32, 
the crop factor for corn is 0.5, and the management 
factor is 0.8 in the upper zone and 0.7 for the other 
zones. 



Results 

Comparisons of the measured and predicted runoff 
characteristics are given in table 15, for the six storms. 
In general, the predicted values exceed the measured 



values with the gieatest error associated with the largest 
storm. The erosion-deposition chatacteristics for these 
six storm events arc given in table 16. The only large 
discrepancy between measuicd and picdicted sediment 
yield is for the June 29 storm when the hydrologic 
characteristics were also most in error. The detailed 
erosion results of each compartment are given in table 
17 for the second storm on June 28, 1965. These results 
illustrate the deposition of sediment in the lower 
compartment of the D stream tube 




Figure 18 -Treynor, Iowa Watershed 2. (82.8 acres) 



Figure 19. -Treynor Watershed W-2 transformed. 
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Table 15,-RunojJ characteristics, Treynor Watershed W-2, 6 storm events- 6-2 '7-65 to 

7-1-65 



Runoff 


Date Ram rail 


Measured 


Predicted 


Volume 


Peak 


Volume 


Pcdk 


6-27 


In 
53 


In 
0.04 
.16 
.68 
.64 
.13 
.05 


Cfs 

20 
72 
152 
157 
51 
21 


hi 

016 
33 
92 
1,51 
.12 
02 


Cfs 
51 
132 
204 
384 
39 
3 


6-28 
6-28 


.62 
1.20 


6-29 


2 42 


7-1 


68 


7-1 


.25 



Table 16. -Erosion-deposition characteristics, Treynor Watershed W-2, 6 storm events- 6-27-65 to 7-1-65 



Date 


Rainfall 


R 


Measured 


Predicted 


Erosion 








amount 




yield 


yield 


Internll 


Rill 




6-27 


In 

0.53 


8 8 


T/ac 
51 


T/ac 


T/ac 


T/ac 


T/ac 


6-28 


.62 


8 5 


1 84 










6-28 


...... 1.20 


31 8 


5 36 


5 74 








6-29 


2.42 


26.4 


3.65 


8 21 


1 01 


1 95 




7-11 


.68 


1.6 


.58 


33 


07 


30 


03 


7-1 


.25 


.5 


.15 


.05 


03 


02 


004 



















Table 1 7. -Erosion results in each compartment for the 2d storm on June 28, 1965 



Compartment 


Interim 




Rill 




Deposition 






Amount 


Area 


Depth 




Al 


Tons 
26.8 


Tons 
38.0 


Percent 
17 7 


In 

0065 


Tons 

o 


Bl 


4.7 


6.7 


17.7 


065 


o 


B2 


4.1 


63.1 


24.6 


.528 


o 


Cl 


39.8 


192.9 


24.6 


.248 


o 


Dl 


11.5 


55.6 


24.6 


.248 





D2 


8.6 


24.8 


20.0 


.170 


43.5 


El 


17.4 


24.8 


20.0 


.064 
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FUTURE PROSPECTS 



As indicated in the discussion of the submodels, the 
hydrology submodel has been tested separately on 
several watersheds and found to produce volumes and 
peaks corresponding reasonably well to the obseived 
data. The coupling of the hydrology and erosional 
submodels lias been tested twice, once with the five 
storms on the Coshocton, Ohio watershed and with the 
six storms on the Treynor, Iowa watershed. The 
combination of all three submodels has been tested once 
with the five storms on the Coshocton, Ohio watershed. 

An adequate test of this model requires complete sets 
of chemical, erosional, and hydrologic data for several 
diverse watersheds. This information is extremely rare 
because it is only lecently that the chemical and physical 
data have been collected simultaneously. The only pub- 
lished data that approach the necessary requiiements arc 
the pesticide work on the Coshocton, Ohio watersheds. 
Pesticide experiments recently completed on theRiesel, 
Texas watersheds and currently in progress on the 
Watkinsville, Georgia watershed may provide some 
additional tests of the model 

The best tests of any model come Irom experiments 
designed for that purpose. Therefore, it is hoped that the 
chemists and engineers at all of our watershed centers 
can use this model as a research tool for designing 
experiments and data collection, The experimental 
results will provide tests for this version of the model, as 
wel! as information for developing better versions. 

These experiments will require an unpredictable 
period of time to complete because of erratic weather. 
In the meantime, the literature is being examined for 
data that might serve to test one or two parts of the 
model. A number of papers have been published on the 
chemicals in runoff and erosion from simulated rainfall 
experiments. Some of these may have sufficient data to 
test the combination of erosion and chemical submodels 
for a single compartment. A number of leaching experi- 
ments might be used to test the leaching aspect of the 
hydrology and chemical submodel combinations. A few 
of die numerous nitrogen experiments might have 
sufficient data to test some aspects of the hydrology and 
nitrate submodel combination. 

It is anticipated that the testing will result in some 
model concepts being revised or replaced. It is hoped 
that this model will serve as a stimulus for conducting 
small experiments to test parts of the model; leaching, 
dispersion, rill development, etc. The development of 
the submodels revealed a number of areas where basic 
research is needed to improve our understanding of 



important processes The following few paragraphs 
highlight some of the most obvious problems. 

The path that the water follows through the soil takes 
on added importance when the chemical aspect is 
considered. The lateral flow of water below the surface 
is extremely difficult to monitor and descube. Interflow 
appears to be the best explanation for the highest 
chemical concentration in runoff occuriing in the second 
and third storms of some expenments. Snow distribu- 
tion and melting are poorly understood, paiticulaily in 
the piediction of runoff and leaching of chemicals. 
Considerable information and a number of evapotrans- 
piration models have been produced recently, but the 
root growth and water extraction pattern in the soil is 
still inadequately understood. 

The dynamic 01 time-dependent aspects of erosion 
take on added importance when chemical adsorption 
and distribution is considered. The detachment, tians- 
poit, and deposition piocesses depend upon rates as well 
as volumes of lainfall and iimoff. The changing rates of 
watei and soil movement will be i effected in time- 
dependent chemical movement. 

Enrichment ratios indicate the need to chaiactcrize 
the sediment movement beyond amount and time. Even 
particle size, while very important, is not sufficient to 
account for the chemical transport properties of ag- 
gregates and organic matter. If the concepts of rill and 
inteinll erosion used in this model are shown to be 
important characteristics of chemical loss by erosion, 
then a great deal of research will be needed to 
adequately describe such processes as the formation and 
configuration of rills. 

In the chemical area the washoff of pesticides from 
vegetation and the leaching of nutrients out of the 
vegetation are some of the most poorly understood 
processes. Management practices, such as banding or 
repeated applications, appear to cause greater effects 
than expected. The immobilization and denitnfication 
of nitrate continue to be grey areas needing research. 
The chroma tographic theory used in this model is being 
tested in the laboratory and in small plots and if found 
to be adequate, considerable work will be needed to 
measure and predict dispersion distribution coefficients. 

The limited target of sloping Corn Belt watersheds 
must be broadened to include all important agricultural 
lands. It is believed that the submodels are adequate for 
most soils. The soils with large cracks are expected to 
cause the most serious problems. Application of the 
model to flat or level lands, such as in parts of the Corn 
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Belt, the delta, and coastal lands, may require alternate 
submodels. The use of large amounts of chemicals in 
these areas makes them a most important challenge. 
Irrigated lands are also flat and receive considerable 
amounts of chemicals. Only parts of the present hy- 



drology and erosion submodels may be needed to model 
the flat and irrigated lands. Salinity has not been 
included because it was felt that there were several 
modeling attempts already underway and thus more 
could be accomplished in these other areas. 
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